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ABSTRACT 

Certain  measurements  of  the  electrical  and  mechanical 
properties  of  transducers  are  of  importance  to  designers. 
These  can  be  determins-d  from  accurate  compls-x  admittance  and 
impedance  measurements  made  throughout  the  frequency  range 
of  interest.  Manual  collection  of  needed  data  is  a 
time-consuming  process  that  Is  prone  to  error.  The 
computerized  system  herein  described  substantially  reduces 
the  time  requirement  and  produces  more  accurate  output  than 
can  be  obtained  utilizing  manual,  methods.  In  addition  10  a 
general  discussion  of  equivalent  circuits  and  the 
instruments  employed,  this  report  includes  samples  of  plots 
and  calculated  parameters  for  piezoelectric  and 
magnetcstrictive  transducers,  the  experimental  comparison 
with  traditional  manual  methods  and  the  program  listing  for 
two  versions  of  tha  system  illowir.g  varying  amounts  of 
operator  options. 
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I.   INTRODUCTION 

Transduction  is  defined  in  an  extremely  general  sense  as 
any  conversion  of  oie  form  of  energy  into  another.  The 
devices  responsible  for  this  conversion,  in  either 
direction,  are  transiicers.  For  this  author,  the  interest 
is  the  conversion  of  slectricity  into  sound  and  vice  versa. 
Commonly,  a  transducer  that  converts  electrical  energy  into 
sound  is  called  a  projector,  loudspeaker,  source  or 
transmitter.  Those  serving  the  reverse  function  are 
commonly  called  hydrophones,  aicrophor.es,  or  receivers. 
Seme  types  of  transducers  are  reversible  and  are  capable  of 
filling  either  role.  This  study  is  restricted  to  the 
analysis  of  reversible  transducers. 

A.   MECHANISMS  OF  TRANSDUCTION 

There  are  two  basic  types  of  mechanisms  commonly  used  in 
reversible  electroac oustic  transducers  designed  for 
underwater  sound,  one  involving  interaction  between 
mechanical  motions  and  forces  and  electrostatic  fields  and 
the  other  between  necaanical  notions  and  forces  and  magnetic 
fields.  Both  types  are  aneaable  to  analysis  jsin.g 
equivalent   electrical  circui-s;    but,    due   to  the   basic 
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differences  in  conventions  used  to  describe  forces  in 
electric  and  in  magnetic  fields,  the  details  are  different. 
Magnetically  coupled  transducers  of  the  moving  coil  type 
have  the  diaphram  attached  to  a  wire  coil  suspended  in  a 
magnetic  field.  It  is  the  interaction  cf  the  current  in  the 
coil  and  the  magnetic  field  that  induces  a  force  on  the  coil 
and  causes  the  diaphram  to  oove  and  radiate  an  acoustic 
wave  [ Ref .  1]-  rhase  classifications  may  be  further 
subdivided  based  on  the  motor  mechanisms  involved  (i.e. 
electrodyna mic,  electrostatic,  magnetic,  aagnetostrict ive, 
and  piezoelectric)  [Bef.  2].  The  classification  systems  for 
transducers  vary  depending  upon  the  analysis  being  done. 
Descriptions  of  variDis  classifications  are  presented  by  the 
previous  references  as  well  as  by  Heaslip  [Ref.  3],  and 
Sherman  [Ref.  4  ]. 

B.   PURPOSE  OF  THIS  STUDY 

All  electroacoust ic  transducers  may  be  considered  as 
being  composed  of  eLactrical  elements  so  that  the  analysis 
of  the  transducer  is  greatly  simplified  if  an  eguivalent 
electrical  circuit  is  used. 

While  the  most  oommon  transducers  currently  in  use  in 
underwater  sound   ars  made  of  piezoelectric   materials,   the 
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aim  of  this  study  is  to  provide  a  generalized, 
semiautomated  means  of  obtaining  meaningful  transducer 
measurements  on  any  type  of  transducer  with  the  operator 
having  mininal  knowledge  of  the  device. 

The  equipment  set-up  utiLizes  basic  devices  (e.g. 
voltmeter,  sweep  oscillator,  Dranetz  impedance  meter, 
amplifier)  interface!  to,  ani  controlled  by,  a  desk-top 
computer.  This  interfacing  provides  more  precise  data 
acquisition  than  can  be  obtained  from  manual  plots 
interpreted  by  an  operator.  The  time  saving  factor  is  a 
major  motivation  for  this  choice. 

Certain  measurements  of  the  electrical  and  mechanical 
properties  of  the  transducer  irs  of  importance  to  designers. 
In  some  cases,  it  is  desireable  to  have  the  resonant 
frequency  fall  within  a  particular  range.  More  likely, 
there  will  be  concern  for  a  high  electromechanical  coupling 
coefficient.  Efficiency  nay  be  the  major  concern  or  perhaps 
it  is  desired  to  have  high  losses  and,  therefore,  a  low 
mechanical  quality  factor  witi  consequent  "flatness"  of 
response  [  Ref .  5  ]. 

Direct  measurement  of  these  values  is  net  always 
possible.    However,  many  of  tha  electrical  component  values 
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needed  to  calculate  lechanical  parameters  can  be  readily 
obtained  from  propacLy  annotatsd  impedance  and  admittance 
plots  for  a  particular  transducer. 

Currently,  there  are  automated  systems  designed  to 
analyze  transducers,  such  as  tis  WQM-12  [Ref.  6 ]•  Few  are 
small  enough  to  be  readily  transported  to  the  transducer  as 
is  this  system, 

C.   FORMAT  CF  THE  REPORT 

The  basic  physical  and  electrical  theory  used  in 
preparation  of  this  report  is  presented  in  Chapter  II. 
While  more  thorough  analyses  of  transducers  and  equivalent 
circuits  may  be  found  in  some  of  the  references,  only  data 
pertinent  to  the  measurements  obtained  by  this  system  are 
presented. 

Chapter  III  expLains  the  interfacing  between  various 
components  of  the  equipment  set-up  and  the  computer.  The 
operator/computer  interaction  is  also  addressed. 
Representative  samples  of  computer  outpu-  are  included. 

Chapter  IV  briefly  describee  the  testing  and  evaluation 
of  the  system  and  program.  Ths  results  and  conclusions  are 
presented  here. 
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The  computer  programs   urittsn    in    BASIC   computer   language 
are   included  in    Appendix      A. 
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II.   THEORY 

Since  much   of  the  physics   and  engineering  involved  in 

discussions    of   transducers    is    common  to    both 
electric-coupled  and   magnetic-coupled   transducers,    this 

report   will   make   a   general    approach   and  point   out 

differences  where  they  occur.    Table  I,      at  the  end  of  this 

chapter,   contains  a  listing  of   symbols  used  in  this  report 

and  refers  to   symbology   used  in   some   of  the  references 
cited. 

A.   GENERAL 

The  approach  usually  taken  in  electroacoustic  analysis 
is  to  use  an  electric  circuit  to  model  the  transducer.  In 
simplest  form  one  may  use  a  two-port  network  with  one  port 
representing  the  electrical  terminals  and  the  other  the 
mechanical  terminals.  Analogs  for  force  and  speed  Nay  be 
voltage  and  current  or  vice  versa  depending  on  the  type  of 
transducer. 

The  canonical  equations  inscribing  the  behavior  of  the 
two  port  transducer  may  be  written  as: 

T'=V  x  *  Tem'U  (2.1) 

and 
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F^me1  +     Zm-°  (2.2) 

where  V  and  I  are  the  voltage  md  current  at  the  electrical 
port  and  F  and  U  are  the  force  and  velocity  at  the 
mechanical  port.  rem  is  a  transduction  coefficient  relating 
the  electromotive  farce  at  tie  electrical  port  to  the 
velocity  in  the  mechanical  net,  and  Tme  is  a  transduction 
coefficient  that  relates  the  force  developed  at  the 
mechanical  port  of  a  two-port  network  to  the  current  in  the 
electrical    mesh   [Ref.    8].  4>  is    the      transformation   factor 

and  may  be  defined  as  4>=Tem/Z0.  The  canonical  equations  for 
a  maqnetically  coupled  transducer  of  the  moving  coil  type 
are: 

F=   -B-l-  I    ♦   Zma  (2.3) 


V-ZQ-  I   ♦    B-1-0  (2.U) 

where  3  is  the  magnetic  field  and  1  is  the  length  of  the 
wire  in  the  coil.  (Mote:  This  choice  is  valid  for 
electrodyna  mic  transducers  (ooving  coil)  as  written  using 
B-l=Tem.  ?lagnetostcictive  transducers  obey  the  same  set  of 
equations  but  Tem#  3-L.) 
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For  electromechanical  coupling,  Figurs  2.1  shows  typical 
circuits  and  transf ocmations  "Rsf.  7]. 


o- 

Y   I 


u 


Ten 


-T, 


en  [ 


\a 


m 


^m 


IS 


MSth  2m 


F    = 


<KI3HSK5 


Fig.  2.1.   Typical  Circuits  and  Transformations 


In  the  vicinity  of  resonance,  which  is  of  concern  here, 
the  Mason  circuit  or  modified  Van  Dyke  circuit  is  valid.  It 
is  also  valid  for  low  frequency  values  which  will  also  be 
needed  for  some  calculations  \  8ef .  8].  This  circuit  is 
shown  in  Figure  2.2  where  Lm  ,  Rm  ,  and  Cm  represent 
mechanical  values.  ZQ  is  the  clamped  or  blocked  or  shunt 
capacitance . 
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Figure  2.2.   (a).  Moiifiei  Van  Dyke  or  HasDn  circuit  (b) . 
Basic  Transducer  Circuit. 

Kinsler,   and  others  [Ref.  1],    define  the   following 
measurable  mechanical  and  electrical  impedances  of  a  system: 


Blocked  electrical  impedaace  (ohms)  =ZQ  =  V/I  !tj_q 


5) 


Free  electrical  impedance  (ohms)=Z  =V/I| 


F  =  0 


(2.6) 


Open    circuit    mechanical   iapedance    (N  •  s/m)  =Zm=F/U| 


(2.7) 
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Short  circuit   mechanical   impedance    (!*•  s/m)  =  Zmm=P/U  lv_0 

(2.8) 

Magnetically  coupled  transducers  may  be  represented  by 
the  circuit  of  Figuca  2.  3,  whsre  $m  is  a  transformation 
factor.      <£m  =   B-l  for    moving   coil    transducsrs   [Bef.    2]. 
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Figure  2.3.      Equivalent  Circuit   for   a   Magnetically  Coupled 
Transducer. 

Important  properties  of  tha  transducer  may  be  found  by 
studying  the  electrical  impedance  of  the  system  as  a 
function    of   driving    frequency.  The    unloaded  driving    point 

electrical    impedance      may    be      determined    from      the   canonical 
equations    (2.1-2. U)     aid   is    given    by 


Z=V/IlF=0=Zo    *    (-Tem^me/Zm) 


(2.9) 
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Similarly,  the  short  circuit  mechanical  driving  point 
impedance  is  given  by 

The  electrical  impedance  is  influenced  by  the  motion  of  the 
coupled  mechanical  system,  as  indicated  by  the  second  term 
on  the  right  hand  siie.  This  term  is  referred  to  as  the 
motional  impedance,  and  defined  by  Hunt  [Bef.  2],   as 

Zmot=  <-TemTme>  <  VZJ  (2.11) 

Similarly,  the  mechanical  driving  point  impedance  is 
influenced  by  the  electrical  load. 

B.   USE  OF  COMPLEX  13?  EDANCE/ADS  ITTANCE  PLOTS 
1-   The  Resonance  Circle 

Graphical  displays  of  inpedance  aid  admittance  data 
can  be  extremely  useful  in  the  analysis  of  a  transducer. 
When  properly  scaled  and  annotated,  quantitative  results  may 
be  obtained  from  the  diagrams.  Generally,  in  the  vioinity 
of  resonance,  the  element  is  regarded  as  having  a  single 
degree  of  freedom  aid  the  locus  of  points  of  the  complex 
impedance  is  a  circls  in  the  nsighborhood  of  the  resonance 
frequency  [  Ref -  7],  [Ref.  3].  Figure  2.4  and  Figure  2.5 
are  typical  resonance  circle  diagrams  where,   in  Figure  2.4, 
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the  electrical  reactance  is  plotted  as  a  function  of 
electrical  resistance  with  frequency  as  a 
parameter  [Bef.  8].  FU  and  f5  are  the  frequencies  of  the 
half  power  points  and  f8  is  the  resonance  frequency.  Dashed 
lines  show  the  same  transducer  inder  load  conditions. 


h 


Fig.    2.4.       Complex    lapedance   Diagram 


(Note:  Reactance      nay        be      entirely        negative      for        an 

electrically  coupled   transducer). 


23 


Fig.  2.5.   Cooplex  Admittance  Diagram 


In  Figure  2.5,  f2  and  f3  equate  to  the  frequencies 
of  the  half  power  points.  FO  is  the  frequency  of  resonance. 
The  dashed  line  would  be  the  same  ideal  transducer  tasted 
under  load  conditions.  (Note:  Susceptanoe  may  be  entirely 
negative  for  a  magnetically  couplad  transducer . ) 

Motional  impedance  components  may  be  obtained  using 
the  following  formulae: 
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Rmot    =RE    "   R« 


(2.12) 


Xraot     =XE    +    ^o 


(2.13) 


with   Zmot       =  Rmot      ♦    JXmot 

The    motional  admittance   componeits   are   given    by 
Gmot      =    GE   -     (1/R0) 


(2.14) 


(2.15) 


Bmot      =    BE  ~    3  u  co 


(2.  16) 


with    Ymot       =  3mot      -    j  Bmot 


(2.17) 


The  frequency  of  maximum  power  output  at  constant  voltage 
input  is  called  the  mechanical  resonance  ind  is  represented 
in  Figure  2.6  by  f8  [Ref.  9  ]- 

In  the  motional  admittance  diagram,  Figure  2.7f  fO 
is  the  frequency  of  mximum  coidnctar.ee  and  resonance.  F2 
and  f3  are  the  half  power  points. 

2.   Com22H§2.i  &§§&&&££§ 

Impedance  measurements  are  made  it  the  electrical 
terminals  of  the  transducer.    The   behavior  of  the  real  and 
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/'  Kmot 


Fig.    2.6.      Motional   Impedance   Diagram 


imaginary  components  of  the  complex  impedance  may  be 
observed  as  the  frequency  spectrum  is  swept  at  a 
sufficiently        slow        rats      so      that  the  system  is 

quasi-stead y-state.  When  the  nechanicai  system  is  clamped 
so  that  it  cannot  move,  the  blocked  electrical  impedance  may 
be  measured  and  used  to  calculate  motional  impedance.  Since 
it  is  difficult  to  adequately  clamp  a  device  over  the  ea-ire 
frequency  spectrum,  values  may  be  obtained  far  above  and 
below    resonance    and    tie   intermediate    values    inferred.         This 
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mot 


Fig.    2.7.      Motional   Admittance   Diagram 


is  especially  important  for  the  reactance  and  susceptance 
curves  whose  clampad  values  vary  more  drastically  with 
frequency  than  do  the  rsal  components  Df  admittance  or 
impedance. 

The  maximum  value  of  notional  impedance  will  eauate 
to  the  mechanical  rasistanca  and  is  the  diameter  of  the 
circle.  This  valua  *ill  be  of  importance  in  calculating  the 
efficiency    of   the  system. 
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Plots  of  various  input  electrical  parameters 
(conductance  (GE)  and  susceptaace  (BE)  or  resistance  (R£) 
and  reactance  (X E )  versus  frequency)  can  provide  all 
of  the  information  necessary  to  obtain  the  electrical  and 
mechanical  quality  factors  far  the  system.  These  values 
indicate  the  sharpness  of  the  resonance. 

a.   Impedance 

The   circuit   of   Figure    2.8    is   useful    for    analysis 
of   the   behavior    in   the    vicinity    of   resonance. 
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Fig.    2.8.      Equivalent    Circuit   Diagram 


A   plot   of      the    input    electrical    impedance      components   versus 
frequency    will    resemble   Figure    2.9. 
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Fig.  2.9.   Resistanca  and  Reactance  Plotted  Versus  Frequency 


b.   Admittanrs 


The  motioaal  admittance  equations  are 


Ymot  =  Gmot  *  3B 


mot 


R  mot  ♦  J  ^ot    B  mot"^  J^ot 


Bmot 


-  X 


^ot  +Xmot 


♦J- 


mot 


Rmot  +  xmot 


(2 .18) 
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Plots  of  conductancs  and  susc=ptance  versus  frequency  for 
Figure  2.8  will  resenble  Figure  2.10.  The  frequencies  where 
the  susceptance  (B)  is  zero  correspond  to  the  electrical 
resonance  and  antiresonance  frequencies. 
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Figure  2.10.   Conductance  and  Sasceptance  Plotted  Versus 
Frequency 

C.   EVALUATION  OF  TRANSDUCER  PARAMETERS 

This  system  is  designed  to   Dbtain  useful  information  on 
devices  used  in  transmitting  iai   receiving  in   water.    As 
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such,  it  is  first  necessary  to  obtain  data  in  an  unloaded 
situation  in  order  t3  be  able  to  separate  the  effects  of 
acoustic  loading.  For  transducers  designed  for  use  in 
water,  air  normally  provides  an  adeguate  medium  for 
evaluation  cf  the  acoustically  unloaded  properties.  Placing 
the  transducer  in  water  leans  the  two-port  network  will  now 
reflect  the  sum  of  the  mechanioal  and  the  load  impedances. 
Then,    according   to  Hunt  [fief.    2], 

ZT  =    Zm+    zl    =    (Rm>Rl>     +    :)<xm*    X,)  (2.19) 

or 

Zm+     Z|    =     (Rm  f    Ri)C    ♦    J-2Q-p)  (2.20) 

where  p  is  a  freguenoy  parameter--  (p=.  5  ( (<J/u0  )-  (u>0Aj)  )  a  nd  Q 

is  the  guality  factor  measur=i  for  the  resonance  (i.e. 
electrical  or  mechanical  depending  on  the  type  cf 
measurements)  . 

Data  obtainable  from  ainittance  neasurements  will 
provide  needed  valies  for  determining  the  motional 
admittance  circle.  The  maximum  value  of  conductance  will 
occur  at  the  freguanoy  of  elactical  resonance.  The  guality 
factor  can  be  determined  from  the  freguerioies  at  whioh  the 
conductance   assumes  half   the   iii:f f erer.oa   of  its   maximum 
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value  and  its  blocked  value  (i.e.  Q  =  f  /(fu  ~fi)«  These 
frequencies  will  also  coincide  with  the  frequencies  (f2  and 
f3)  of  the  local  maximum  and  minimum  susceptance  values  near 
resonance.  Blocked  capacitance  may  be  calculated  from 
measurements  of  susceptance  far  above  and/or  below 
resonance.  For  an  electricaLly  couplel  transducer  with 
small  dielectric  losses,  valid  measurements  of  C0  may  be 
made  very  far  below  resonance.  For  magnetic  coupling, 
measurements  must  be  taken  far  above  and  far  below  resonance 
and  the  values  between  interpslated  linearly  between  these 
two. 

Impedance  data  will  provide  the  frequency  of 
antiresonance.  The  frequencies  of  the  half  values  for  the 
difference  between  tie  maximum  resistance  and  the  blocked 
resistance  will  correspond  to  tis  maximum  and  minimum  values 
of  reactance  and  allow  for  calculation  of  a  mechanical 
quality  factor.  Blocked  capacitance  or  inductance 
(depending  on  the  typs  of  transducer)  may  be  calculated  from 
reactance  data  obtained  far  abov=  and/or  bslow  resonanci. 

The  following  paragraphs  briefly  describe  the  definition 
and  calculation  of  various  constants. 
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1-      Quality   Factor    =   2 

This  factor   indicates   the   sharpness    of   resonance. 


frequency   of    max.    S   or    R               f (resonance) 
Q   ■ ■» ■    

freq.    difference    of   mm.  f  (upper)  -f  (lower) 

and  max.    B   or   X 

(2.21) 


2.      Electromechaai  cal   coupling   coefficient   f  5.£(§If) 

This        coefficient        indicates        how        tightly        the 
electrical    and    mechanical  meshes   are   coupled   [Ref.    2]. 

K2  (ef  f)    =    1- (f3/f0)  2  for    electromagnetic   coupling 

(2.22) 
and 

K2(eff)    =    1-(f3/f8)2  for    electrostatic   coupling 

(2 .23) 

where      f8      is   used      to      indicate      the  frequency      of      maximum 

resistance      and      f3      indicates        the  frequency      of      maximum 
conductance. 

3-      Static    Coupling.  Coefficient    =  K£ 

This     term    is      the    ratio      of  the   stored      mechanical 
energy   to      the    total    energy      stored    in      the    device  [Ref.     1  ]v 
[Ref.    5]. 
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cl 

K2- (2.  2ft) 

Cn  ♦  C, 


ft.   Shunt  Capacitance  *  C0 

This  is  ths  blocked  capacitance  of  the  system. 
There  is  no  simple  nethod  of  measuring  this  quantity  (with 
the  possible  exception  of  instruments  designed  for  a 
specific  frequency).  fleasur sments  must  be  made  at  a 
frequency  far  removed  from  frequencies  of  mechanical 
resonance  in  order  to  ninimize  the  effects  of 
electromechanical  coupling.  If  the  capacitance  ratio 
(C;/C0)  is  greater  than  fifty,  the  inaccuracy  should  be 
minimal  and  may  be  considered  negligible. 


1       3 
C0  = = where  <o   <<  o>0  (2.25) 

o>XG      oj 


5.   Clammed  Resist  a  nee  and  R3act1ar.ee  =  Ec  ^£.2  £q 

R0  and  X0  may  be  closely  approximated  by  obtaining 
data  far  below  and  above  resonance  and  interpolating  between 
since  there  is  a  positive  frequency  dependence.  However, 
for  many  practical  purposes,  (and  for  elsctricaily  coupled 
transducers),  data  may  be  obtained  far  below  resonance  which 
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provides  a  reasonable  approxiaation  to  the  actual  value 
throughout  the  frequency  range  of  interest.  These  values 
are  assumed  to  represent  those  which  would  be  measured  if 
the  mechanical  parts  of  the  system  were  blocked  or  clamped 
to  prevent  motion. 

6-   Blocked  Inductance  =  Lq 

This  is  the  inductance  of  the  blocked  system 
(usually  with  magnetic  coupling)  and  is  obtainable  from  data 
collected  far  below  and  above  rssonance. 


L     = 


X(interpola  ted) 


<o. 


(2.26) 


7.   Capacitance  Ratio  =  ]_Rl 

This  is  the  ratio  of  the  blocked  capacitance  (CQ)  to 
the  motional  capacitance  (Ct). 


»R«  = 


PI 
"I" 


":<2 


where  K2  is  the 
static  coupling 
coefficient 

(2.27) 
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8-   ZL2^£§  of  Merit 

The  figure  :>f  marit  is  little  nantioned  in  the 
references  probably  because  it  is  an  arbitrary  formula 
depending  upon  a  designer's  criteria.  Miller  [Hef.  5  ]f 
mentions  K2-Qm  as  ona  possibility  and  utilization  of  the 
coupling  coefficient  as  another.  No  specific  calculation  of 
this  value  is  done  ia  this  rapoct,  leaving  it  to  the  user  to 
insert  as  desired. 
9.   Pi  a  meters 

Impedance  or  admittance  circles  ace  transformsd  to 
motional  diagrams  by  subtracting  blocked  values  or  blocked 
values  multiplied  by  angular  freguency  as  discussed  in 
section  B.  1  above.  The  diameters  of  tha  resulting  circles 
will  provide  needed  values  for  efficiency  calculations.  D^ 
and  D^  signify  the  diametars  of  the  motional  circle  for  the 
transducer  in  water  and  in  air,  respectively. 

D.   EFFICIENCY 

The  overall  efficiency  of  a  transducer  is  the  ratio  of 
the  power  delivered  to  an  external  load  connected  at  the 
output  terminals  to  the  total  pDwer  at  tha  input. 

At  resonance,  this  can  be  easily  calculated  from  known 
values.  For  an  eLectrically  coupled  transducer,  the 
efficiency  is  given  by: 
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D,j.(D4  -Dw) 
=    -*A-  Si-  (2.28) 


S    DA 


where   D4    and  D^   are   the  diameters    of   the   motional   admittance 
circles   and  where  3    is    the   condictance   at   resonance. 

Similarly,  the        efficiency      at        resonance        for        a 

magnetically  coupled  transducer    is   given   by: 

E   a    •JRlJL^.VL.  (2.29) 

R-DA 

where  these  are  diameters  of  the  motional  impedance  circles 
and  R  is  the  resistance  at  resonance  [  Ref  •  10]. 

The  mechanical  power  utilization  factor  is  par-  of  the 
overall  efficiency.  It  is  given  by  R//(R/*Rm)  which  oan  be 
expressed  in  terms  of  loaded  and  unleaded  values  of  Q  or  by 
the  diameters  of  circLes. 


-i =  -— d £-.  =  -Ji. IK...  (2.30) 


Rm  *  R/  °4 


(Note:   All   values  should  be   taken  from   either  admittance 
data  or  impedance  data  but  should  not  be  sized.) 
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Once  the   efficiency   at   resonance   is   determined, 
utilization  of  the  impedance  data   of  a  magnetically  coupled 
transducer  allows  for  calculation  of  the  lagging  phase  angle 
<(3). 


P»  0.5-  ArcCos  ((Rres  -R0|  /Dw) 


(2.31) 


The   frequency   parameter   for   maximum   efficiency   may   be 
calculated  as 


p  =  (DwSin(2'P)  >  /(4-R0.Q 


w' 


(2.32) 


and  from  the  ear  liar  relationship  between  p  and  the 
frequency  ratios,  the  frequency  for  maximum  efficiency  for  a 
magnetically  coupled  transducer  may  be  calculated. 


f  =  f 0  (2p  ♦  (p2*!)1/2) 


(2.33) 


The  frequency  of  maximum  efficiency  may  not  coincide 
with  either  the  resonant  frequency  or  the  frequency  cf 
maximum  admittance.  For  the  magne+ostrictive  transducer, 
this  frequency  of  maximum  efficiency  is  usually  greater  than 
that  at  resonance. 
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A  development  by  Camp  *  Bef •  8]f  shows  that  a 
magnetcstr ictive  transducer  may  produce  aiDce  acoustic  power 
at  the  resonant  frajuency  even  though  efficiency  may  be 
less.  Normally,  the  fraquancy  for  optimum  acoustic  power 
will  fall  between  ths  mechanical  resonance  frequency  and  the 
frequency  of  optimum  efficisnoy.  This  is  due  to  eddy 
current  losses  and  magnetic  hystarasis. 

Since  transduction  losses  a:=  lowr  maximum  power  output 
for  piezoelectric  iavicas  is  limited  mainly  by  voltage 
limits.  The  optimum  operating  fcsquency  is  accepted  as  the 
frequency  of  maximum  admittance  [  Ref  •  8]. 

If  one  assumes  that  the  transducer  has  a  high  Q  (sharp 
resonance)  and  that  the  blockad  resistance  is  constant 
throughout  the  fraquancy  ranga  for  the  impedance  circle, 
then  the  potential  efficiency  may  be  expressed  as 


Pot.  Eff.  =  -iasu      -IswlL.. (2#34) 


The  potential  efficiancy  reprasar.ts  the  maximum  efficiancy 
available  for  the  most  favorabla  loading  conditions.  This 
equation   deletes  tha   requiramant  for   data   from  a   loaded 
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condition  but  requires  the  air  motional  impedance  loop  to 
very  closely  approximate  a  circle  [ Ref .  9].  "High  Q "  as 
used  here   implies     values  greater    than   thirty.  If   the    air 

motional  impedance  loop  is  not  approximately  circular  and/or 
the  Q   is   low,    then,    from  refersace    2, 


(1    ►  §*  COS«0  »  V*    -     (1-§»SIN2p)V* 


(1     ♦    }£   C0S2  £  |   V*     -f      (1-  J£    SIN2  j3    )  V2 


(2.35) 


E.       SUMMARY 

There    is      helpful    informatiDn    to      be   obtained      from    both 
impedance      and      admittance      data.  Some      transducers      more 

readily  reveal  properties  through  one  type  of  diagram  than 
another.  Normally,         admittance      data      is      collected      for 

electrically  coupled  systems,  while  impedance  data  is  better 
for  studying  magnetically  couplsi  systems  [Ref.  2].  Due  to 
the  ease  of  obtaining  data  with  the  system  proposed  here, 
both   sets   of  data   ani    diagrams    nay    be    rapilly    obtained. 
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III.       ADAPTATION    TO    A    22MPOTERIZED    SYSTEM 

In-phase  and  out-of-phase  components  of  the  admittance 
or  impedance  are  needed  in  tha  analysis  of  a  transducer.  An 
initial  attempt  was  nade  to  ass  a  lock-in  analyzer,  but 
ultimately    a  Dranetz   impedance   leter    was   used. 

A.       EQUIPMENT    SET-OP 

Data  needed  for  the  calculation  of  properties  of 
transducers  may  be  obtained  usiig  a  lock-in  analyzer.  The 
magnitude  and  phase  or  in-ptiase  and  quadrature  components 
will  indicate  resonances  aid  provide  the  necessary 
information. 

The  initial  set-up  utilized  the  Princeton  Applied 
Research      Company      Model      5204         Lock-in      Analyzer.  Phase 

information  would  have  allowed  for  desired  calculations.  In 
experimentation,  the  circuit  required  for  obtaining  valid 
phase  information  became  too  coaplex  to  be  practical.  The 
lock-in  analyzer  is  designed  to  operate  in  the  1  Hz  to  100 
KHz  frequency  range  although  it  was  determined  to  be 
accurate  at  higher  frequencies.  At  low  frequencies,  there 
is   a      low   source      resistance   and    the      signal    to      noise    ratio 
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became  too  small  to  provide  a  reference.  It  was  concluded 
that  complex  impedance  measursnents  would  be  another  means 
of  getting    the   necessary  data. 

As  briefly  mentioned  in  the  introduction,  there  are  many 
methods  for  obtaining  impedanca  information.  Since  this 
au+hor  desired  a  portable  systsm  usabla  for  underwater 
transducers,   a    bridgs    network   sssmed   reasonable. 

The  set-up  used  to  design  tiis  system  is  shown  in  Figure 
3.1.  The  Hewlett- Pa o kard  HP-35  computer  directs  the 
peripherals  for  the  data  talcing.  The  printer  and  plotter 
are  non-essential  to  the  systsm  and  were  used  to  provide 
larger  format  graphics  and  printout  than  are  available  with 
the   built    in  thermal   printer   provided    with    the    HP-85. 

1.      Hewlett-Packai:  d   HP-85    Computer 

The  HP-85  is  an  eight  bit  microcomputer  that 
utilizes  BASIC  computsr  language.  It  has  a  127  millimeter 
diagonal  black  and  white  slectromagnetic-dsf lection  CRT.  A 
32  character  per  line  thermal  printer/plotter  is  part  of  the 
unit.  The  computer  has  16K  bytes  of  read/write  memory  of 
which      14,579   are      available      to      the   operator.  For      this 

program  a  memory  module  is  nssded  to  expand  the  memory  to 
32K    bytes.      Programs   or  data    ma/    be    stored    on,    or    read   from, 
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magnetic      tape    cartridges.  To      interface    with      peripheral 

equipment,    an    I/O  ROM      and   an    iiterface      card   were     added  to 
provide        HP-IB        (IESE-438)  instrumentation  interface 

capabilities. 

2-      Syn thesiz5r/Fii nctior.   £eaerator 

The  Hewlett-Packard  Modsl  3325A  Synthesizer/Function 
Generator  can  produce  five  different  continuous  waveforms. 
For  this  application  a  sine  wave  was  desired  for  which  this 
instrument  has  a  frequency  range  of  1  microhertz  to  20 
megahertz.  Frequency  may  be  specified  with  up  to  eleven 
digits  of  resolution.  Output  amplitude  is  1  millivolt  to  10 
volts  peak-to- peak  into  a  fifty  ohm  load.  This  model  is 
fully  programmable  ttxrough  the  rear  panel  Hewlett-Packard 
Interface  Bus    (HP-IB)  . 

3*      Dri.Ugtz    Comp_lex   Impedan;a--Admittance   Meter 

The  Model  10DC  Complex  Impedance-- Admittance  Meter 
(CIAM-100C)  is  a  transistorizsd  instrument  that  measures 
vector  impedance,  vactor  admittance,  vector  amplitude  and 
phase.  It  covers  a  frequency  range  froa  100  Hz  to  233  KHz 
in  three  ranges  of  3.1  to  2  KHz,  1  to  20  KHz  and  1  0  to 
200KHz.  For    measuring      impedance      a      constant    current      is 

produced    in      the   unknown   impedance.  The    amplitude      of    the 
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current  is  set  by  tha  range  switch.  The  voltage  developed 
across  the  unknown  impedance  is  amplified  and  fed  to  two 
resolvers.  For  measuring  admittance,  a  constant  voltage  is 
applied  to  the  unknown.  The  current  through  the  unknown  is 
measured  by  means  of  voltage  developed  across  a  small 
range-selected  resistor  placed  in  series  with  the  unknown. 
The  resolver  outputs  are  proportional  to  conductance  and 
susceptance  of  the  n known.  The  claimed  accuracy  of  the 
instrument  is  plus  or  minus  two  percent  full  scale 
amplitude.  The  Dranetz  meter  requires  an  input  signal 
voltage  of  3  Vrms.  In  taking  measurements  on  an  unknown 
transducer,  initially  the  frequency  of  interest  may  not  be 
known  as  well  as  the  appropriate  scale  settings  for  the 
meter.  Measurements  made  on  the  Dranetz  meter  using  a 
variable  resistance  box  and  voltmeters  indicate  linearity 
throughout  the  range  settings.  However,  it  is  necessary  to 
insure  measurements  io  not  cause  either  meter  range  to  be 
exceeded.  The  limit  indicator  Lights  are  activated  when  the 
unknown  impedance  or  admittance  is  substantially  greater 
than  the  full  scale  setting.  ihen  this  occurs,  an  internal 
relay  shunts  the  meters  and  trie  DC  output  terminals  to  zero. 
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U.   Data  Acquisition/Control  Onit 

The  Hewlatt-Packard  Model  3U97A  Data 
Acquis it ion /Control  Unit  is  used  here  as  a  voltmeter  and 
scanner  to  measure  tie  DC  output  voltage  proportional  to  the 
deflection  of  the  meters  of  the  Dranetz.  The  instrument 
measures  and  displays  voltage  values  to  five  and  one  half 
digits.  It  is  internally  triggered  by  software  command 
during  the  data  talcing  but  takes  the  next  reading  as  soon  as 
the  first  is  completed  for  obtaining  averages. 

5«   Power  Amplifier 

This  instrument  is  not  sssential  to  the  operation  of 
the  system.  The  Dranetz  metsr  is  designed  to  accept  direct 
ac  voltage  input;  however,  the  particular  meter  used 
exhibited  more  stability  when  tae  amplifier  was  included  in 
the  set-up.  The  ampLifier  had  unity  gain  and  was  only  used 
to  lower  the  effective  output  impedance  of  the  frequency 
synthesizer  from  fifty  ohms  to  two  milliohms. 

6.   General  Discussion 

This  equipment  set-up  is  suitable  for  an  initial 
system.  A  system  designed  for  regular  and  rcurine  use 
should  include  a  different  imp=dance  measuring  system  than 
the  Dranetz.    Because  of  the  design  feature  of  the  Dranetz, 
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a  great  deal  of  operator  interface  is  required  to  change  the 
various  different  scale  and  raage  factors  which  cannot  be 
set  directly  by  the  computer.  It  is  too  easy  for  an 
operator  to  miss  changing  oae  setting  or  to  overlap 
frequency  range  scales  in  the  spectrum  sweep  and,  thereby, 
obtain  erroneous  data.  IIsd,  the  scale  factors  must  be 
manually  entered  prior  to  the  computer  doing  any  cal- 
culations. The  Dranetz  also  hai  annoying  features,  such  as 
drift  of  the  zero,  excessive  noise  on  low  impedance  ranges, 
and  DC  output  at  full  scale  deflection  of  0.8  -  0.95  volts. 

B.   DISCUSSION  OF  THE  PR03EAM 

The  HP-85  is  the  controlling  unit  directing  the  other 
instruments  to  take  or  send  data  as  needed.  Memory  storage 
is  handled  by  this  microcomputer. 

The  goal  was  to  lesign  a  program  than  would  minimize  the 
time  necessary  to  coLlect  pertinent  data  while  maximizing 
accuracy.  The  program  was  to  accommodate  different  media, 
allow  for  Collection  of  admittance  or  impedance  data  or 
both,  and  work  for  transducers  with  either  electric  or 
magnetic  coupling,  while  being  easy  to  operate.  This  was 
accomplished;  but  due  to  the  large  differences  in  properties 
of  different  transducers,  a  second  program  was  written  tc  be 
used  by  more  experienced  operators. 
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One  program  is  dasigned  so  that  a  person  with  minimal 
knowledge  of  the  transducer  or  measurement  procedures  may 
run  it  and  obtain  meaningful  information.  It  is  suggested 
that  the  operator  bs  familiar  with  the  operation  of  the 
interfacing  instruments  (especially  the  Dranetz  meter).  An 
equivalent  program  is  designed  so  that  the  operator  may 
select  the  bandwidth  for  data  collection  and  may  specify  the 
frequency  fcr  measursment  of  tas  blocked  values  (done  as  a 
percentage    of   the  resonant    frequency) . 

The  program   is      broken    into    three   sections.  The   first 

gives  the  operator  an  overview  of  the  conductance  and 
resistance  of  the  transducer  under  test.  Fart  two  is  the 
data  collection  for  desired  measurements  with  options  for 
plots  and  lists  of  the  data.  The  final  section  performs  the 
calculations  for  motional  lata,  efficiencies,  and  other 
desired   outputs.  The  following      subdivisions   address      the 

program     in   detail.  Figure    3.2      is      a   flow      chart    of      the 

program.  The  programs  written  in  BASIC  are  included  as 
Appendix    A. 

1-      Xk®   Search 

This  first      portion    has    been      included   to      allow    for 
visual   and    graphic   display    of   tie    response   of    the   transducer 
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over  a  selected  frequency  ranga.  The  operator  is  required 
to  insert  information  on  the  type  of  coupling  (magnetic  or 
electric)  of  the  transducer  and  the  medium  (air  or  water)  in 
which  measurements  are  to  be  taken.  The  operator  is 
instructed  to  set  the  Dranetz  for  the  measurement  of 
admittance  data  and  to  specify  the  range  to  be  covered  in 
the  frequency  sweep.  The  operator  also  enters  the  desired 
voltage  to  be  sent  to  the  Dranetz.  This  manual  voltage 
input  is  necessary  since  at  higher  frequencies  lower  rms 
voltages  are  necessary  to  keep  the  ac  voltage  meter  on  the 
Dranetz  at  3.0  V. 

Through  the  Data  Acquisition  Unit,  the  relative 
amplitude  of  the  real  part  of  :he  admittance  is  samplad  at 
300  equally  spaced  points  throughout  the  selected  frequency 
range.  k  graph  is  displayed  on  the  CRT  and  then  copied  on 
the  thermal  printer  of  the  HP-95  for  later  reference.  In  the 
event  that  some  paranater  needs  to  be  changed,  the  operator 
is  given  the  opportunity  to  rerm  this  portion. 

The  operator  is  instructed  to  specify  a  relative 
amplitude  (from  0  to  1.2)  that  is  less  than  the  maximum  of 
the  admittance  peak  d  f  interest.  The  computer  print.5  the 
relative   a  apiitudes    greater   than   this   value    and   the 
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associated  frequencies.  This  allows  the  operator  to  select 
the  resonance  of  interest  and  provide  the  center  frequency 
and  halfwidth  to  be  considarad  in  finding  the  resonant 
frequency  and  quality  factor  Q|  .  A  comparison  is  done  to 
find  the  absolute  maximum  valua  for  amplitude  and  the  half 
power  points  for  the  determination  of  2.  The  resultant 
factors  are  printed  and  the  option  to  rerun  is  again 
extended. 

This  same  subroutine  is  again  used  to  obtain  similar 
data  for  the  real  part  of  the  inoedance.  The  subroutine  is 
an  adaptation  of  a  part  of  a  procedure  designed  by 
Conte  [Ref .  13 ]. 

Although  two  repeats  of  the  subroutine  take  about  10 
minutes,  the  operator  is  given  sufficient  information  on  the 
impedance  and  the  admittance  responses  of  the  transducer  to 
decide  which  type  of  measurements  will  provide  the  most 
accurate  information. 
2«   Da* 3  Collection 

The  computer  has  the  frequency  synthesizer  send  a 
frequency  far  below  resonance.  The  meters  cf  the  Dranetz 
are  zeroed  by  the  operator  and  set  for  full  scale 
deflection.     An   average   of   ten   measurements   for   each 
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54 


I 


PLOT 


Yes 

LIST?  ^~K 


CALCULATE 

MOTICNAL 
DATA 


PRINT 


Figure  3.  2   (Continued) 


55 


Yes 


Figure  3.  2   (Continuod) 


56 


terminal  is  placed  in  memory  to  be  used  as  normalization 
factors  for  all  later  DC  outpit  voltage  measurements. 

While  still  at  this  frequency  far  below  resonance, 
an  average  of  ten  measurements  will  be  made,  normalized  and 
scaled  (using  operator  input  for  the  scale  factors  from  the 
Dranetz  instrument)  and  stored  as  blocked  values  of 
conductance,  susceptaace,  resistance  and  reactance.  These 
will  be  used  later  La  computations  after  being  assigned  to 
different  variables  depending  on  the  medium.  This  one 
measurement  is  all  that  is  necessary  for  the  electrically 
coupled  transducers  with  small  dielectric  losses.  Foe  the 
magnetically  coupled  transducers,  the  values  are  taken  far 
below  resonance  and  stored.  Tie  program  Loops  back  through 
this  section  after  sending  a  message  to  the  frequency 
synthesizer  to  go  an  equal  frequency  distance  above 
resonance.  An  average  of  the  two  values  is  stored  for  later 
use.  This  averaging  is  necessary  since  the  eddy  current  and 
hysteresis  losses  cause  a  steady  increase  in  the  blocked 
impedance  values. 

The  decision  of  tie  operator  en  the  type  of  data  to 
be  collected  is  entered.  The  scale  factor  from  the  Dranetz 
instrument  is  entered  and  the   computer  goes  to  a  subroutine 
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for  data  collection.  The  computer  draws  the  appropriate 
center  frequency  from  the  data  collected  during  the  search 
and  bases  the  bandwidth  for  the  sweep  on  the  quality  factor 
measured  previously.  Fifty  frsguencies  and  the  associated 
real  and  imaginary  lata  are  maasured  and  stored  from  ten 
times  the  half  power  bandwidth  below  center  frequency  to  two 
times  the  bandwidth  below  resonance.  Two  hundred 
measurements  are  taken  in  the  vicinity  of  resonance  and 
another  fifty  in  a  region  from  two  to  ten  times  the 
bandwidth  above  resonance.  (Note:  This  program  reduces  the 
range  of  interest  to  plus  and  minus  five  times  the  bandwidth 
around  resonance  when  tests  are  done  under  load  (water) 
conditions.  The  seoond  progran  takes  data  at  three  hundred 
equally  spaced  frequencies  within  the  bandwidth  specified  by 
the  operator.)  Once  the  data  aave  been  collected,  another 
subroutine  is  entered  to  find  the  minimum  and  maximum  values 
for  the  two  sets  of  iata  and  their  associated  frequencies. 
This  is  accomplished  through  a  comparison  to  a  minimum  (or 
maximum)  value  until  a  smaller  (larger)  one  is  found.  A 
reassignment  is  done  and  the  comparison  continues  throughout 
the  data.  The  maximum  value  and  the  two  adjacent  values  are 
then  fitted  to  a  quadratic  equation  which  is  used  to  produce 


53 


a  better  estimate  of  the  maximum  and  its  associated 
frequency.  In  the  event  of  an  error  in  data  collection, 
there  is  an  opportunity  to  retake  the  data  before  proceeding 
in  the  program. 

A  subroutine  provides  far  a  listing  of  the  collected 
data  (100  points)  in  the  vicinity  of  resonance.  All  of  the 
data  could  be  printed;  however,  the  data  near  resonance 
provides  the  important  information  and  printing  is  an  option 
open  to  the  operator. 

Plots  of  tha  data  may  bs  made.  The  operator  may 
chose  to  make  plots  of  real  or  imaginary  amplitudes  versus 
frequency  or  the  plots  of  real  versus  imaginary  data  over 
frequency,  or  all,  or  none.  Each  plot  is  drawn  using  a 
subroutine  and  labelad  with  pertinent  scaling  and  captions. 
At  the  end  of  this  chapter,  Figures  3.3  through  3.14  are 
examples  of  plots  foe  a  piezoelectric  transducer.  (This 
example  is  for  a  Type  100-635  3-010,  Serial  No.  1213, 
designed  for  use  i.i  an  active  sonobuoy.)  Figures  3.15 
through  3.2  2  are  examples  of  plots  for  a  tnagnetostrictive 
transducer.  (This  sscampla  is  i  Type  CMT-10255,  Model  QGB, 
Serial  No.  318.) 
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3«   Data  Manipulat  ion 

In  the   calculation  of  notional   data,   the  type  of 
coupling  of  the  transducer  becones  important. 

a-   Magnetic  Coupling 

For  a  nagnetically  coupled  transducer,  the 
motional  impedance  data  is  calculated  by  subtracting  the 
blocked  resistance  from  each  resistance  measurement. 
Blocked  inductance  at  resonant  is  calculated  using  the 
averaged  blocked  reactance  value.  The  angular  frequency 
multiplied  by  the  blocked  inductance  is  then  subtracted  from 
each  value  cf  reactance.  In  the  event  admittance  data  has 
been  collected,  the  motional  conductance  is  obtained  by 
subtracting  the  blocked  conductance,  while  motional 
susceptance  may  be  obtained  by  subtracting  the  blocked 
capacitance  multiplied  by  the  angular  frequency  from  each 
value.  A  list  of  notional  data  may  be  obtained.  k  plot  may 
be  elected  and  dra*i  and  labeled  with  axis  scaling  done 
using  minimum  and  maximum  values  of  the  calculated  data. 

Computations  are  dcie  using  the  collected  data. 
The  dynamic  electromechanical  coupling  coefficient  is 
calculated  from  the  resonant  frequencf.es.  k  new  quality 
factor  is  calculated  based  on  tie  frequencies  of  maximum  and 
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minimum  amplitudes  of  the  imaginary  part  of  the  data 
collected.  Values  of  some  blocked  measurements,  quality 
factors,  the  coupling  coefficient,  and  resonances  are  output 
to  the  printer. 

b.   Electric  Coupling 

For  an  electrically  coupled  transducer, 
dielectric  losses  ar=  small  and  it  is  acceptible  to  use  the 
values  measured  far  below  resonance  as  the  blocked 
measurements.  Motional  resistance  is  obtained  by 
subtracting  the  blocked  resistance  from  each  measurement. 
The  addition  of  tha  angular  frequency  multiplied  by  the 
calculated  blocked  iiiuctance  t3  the  reactance  measurements 
will  provide  the  motional  reactance.  The  lata  may  be  listed 
and  a  plot  male.  An  example  of  the  motional  impedance  plot 
can  be  seen  in  Figure  3.9. 

Motional  admittance  data  may  be  obtained  by 
subtracting  the  reciprocal  of  the  blocksd  resistance  from 
each  conductance  valu=.  The  angular  frequency  multiplisd  by 
the  blocked  capacitance  is  subtracted  from  each  susceptance 
measurement.  Figure  3.10  shows  a  motional  admittance  plot 
for  a  piezoelectric  transducer  In  air.  Calculations  similar 
to  those  for  magnetic  coupling  are  done  and  desired  values 
printed. 
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After  measurements  are  taken  with  the  transducer 
under  load  (in  water)  calculations  are  done  to  find  the 
efficiency  and  optimal  operating  frequency  (for  magnetic 
coupling).  This  information  is  output  to  the  printer.  An 
opportunity  to  repeat  data  collection  in  either  medium  is 
afforded  before  lata  on  air  and  water  temperatures  and 
transducer  identification  numbers  are  requasted  prior  to  the 
end  of  the  program.  A  sample  oitput  (for  the 
magnetostr ictive  transducer  Typa  CMT-10255,  Serial  No.  318, 
Model  QGB)  is  included  as  Figura  3.23. 
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VALUES  MEASURED  IN  AIR  FROM  ADMITTANCE  DATA 

BLOCKED  INDUCTANCE  =  8.86E-4  HENRIES  AT  23956  HZ 

ELECTRICAL  QUALITY  FACTOR  =  189.73 

MECHANICAL  QUALITY  FACTOR  =  29.426 

ELECTRICAL  RESONANCE  =  23956.0379323  HZ 

MECHANICAL  RESONANCE  =  23920  HZ 

DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIENT  =  .003 

BLOCKED  RESISTANCE  =  33.724  OHMS 


VALUES  FOR  MEASUREMENTS  IN  WATER  FROM  ADMITTANCE  DATA 

BLOCKED  INDUCTANCE  =  9.01E-4  HENRIES  AT  23932.48  HZ 

ELECTRICAL  QUALITY  FACTOR  =  18.513 

MECHANICAL  QUALITY  FACTOR  =  19.688 

ELECTRICAL  RESONANCE  =  23932.4882297  HZ 

MECHANICAL  RESONANCE  =  23892  HZ 

DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIENT  =  .00338 

POTENTIAL  EFFICIENCY  =  .20398 

FREQUENCY  OF  OPTIMUM  EFFICIENCY  -  24246.61  HZ 

MECHANICAL  POWER  UTILIZATION  FACTOR  =  .14645 


Figure  3.23  Sample  Computer  Output 
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IV.   TEST  AND  EVALUATION 

Although  many  of  the  calculations  mads  by  the  computer 
are  basic,  it  was  desireable  to  have  a  set  af  data  collected 
manually  to  ensure  accuracy  of  the  semiautomated  system. 

A.       TESTING 

1.  Method 

Data  were  collected  for  four  different  transducers, 
all  of  which  were  later  used  to  test  the  computerized 
system.  A  similar  equipment  set-up  was  used  with  the 
computer  replaced  by  an  operator  and  an  X-Y  plotter,  and  a 
Wavetek.  pulse/f unction  generator  in  place  of  the 
synthesizer /function  generator. 

The  frequency  spectrum  was  swept  manually  with  plots 
made  of  the  real  and  imaginary  components  versus  frequency. 
From  these  plots,  amplitudes  and  frequency  values  were 
obtained  and  calculations  similar  to  those  done  by  the 
computer  performed. 

2.  Comparisons 

For  the  magnat ostr ictiva  transducer  used  here,  the 
data  output  of  the  computer  corresponds  very  closely  to  that 
collected   manually.         The    blocfcid    inductance    was    found   to    be 
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the  same  calculated  by  eithsr  method  to  five  digits. 
Resonance  frequency  variations  up  to.  ten  hertz  may  be 
attributed  to  temperature  variations;  and  these  frequency 
variations  will  affect  the  sapling  coefficient.  The 
efficiencies  compared  to  within  five  percent. 

Comparisons  of  the  piezoelectric  data  were  as 
favorable  as  that  for  the  magnet ostrictive  transducer.  The 
blocked  capacitances  varied  by  up  to  six  percent  and  there 
were  small  resonance  frequency  differences  in  some  cases. 

B.   EVALUATION 

The  resonance  frequency  variations  noted  may  be 
attributed  trainly  to  temperature  changes  from  one  experiment 
to  another.  Larger  deviations  noted  in  a  few  cases  are  due 
to  inaccuracies  in  reading  graphs  or  in  determining  a  true 
maximum  for  a  very  Low  3  devioe  as  is  the  case  for  many 
transducers  when  plaoed  under  load.  For  the  examples  used 
here,  the  computer  data  are  considered  more  accurate  since 
it  measures  and  compares  data  to  five  digits,  far  better 
than  the  average  eye  can  do  with  unsealed  data.  In  the 
computerized  system,  the  data  collected  in  the  vicinity  of 
resonance  are  compared  to  find  the  absolute  maximum  value 
and  the  associated  frequency.    This   point  and  two  adjacent 
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points  are  then  used  for  a  parabolic  fit  to  find  the  precise 
frequency  for  maximum  componant  amplitude.  This  provides 
accuracy  in  the  resonance  fraquancy  far  in  excess  of  what 
may  be  achieved  manually  by  an  sperator. 

A  *  delay  time'  to  allow  stabiization  of  the  system 
(dominated  by  the  response  time  of  the  Dranetz  meter)  was 
determined  to  be  oaa  hundred  milliseconds.  In  manually 
sweeping  the  frequancy  spectrum,  stabiization  time  is 
difficult  to  achieve.  This  may  result  in  the  appearance  of 
spikes  in  the  plots  which  incraase  the  difficulty  of  manual 
evaluation  of  the  collected  lata.  Once  this  stabiization 
time  was  determined,  all  wait  times  in  tha  computer  program 
were  ensured  to  be  greater  than  this  value  with  the  longest 
waiting  period  duriig  deteciina tion  of  the  shunt  and 
normalization  values.  The  seoonfl  longest  wait  period  occurs 
during  the  data  collaction.  Hare,  a  period  of  two  hundred 
fifty  milliseconds  has  been  used  in  the  vicinity  of 
resonance  tominimiza  the  chanca  of  error.  The  added  minute 
for  data  collection  that  this  additional  wait  time  requires 
is  of  minor  inconvenience  compared  to  tha  improvement  in 
accuracy. 
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Transducers  may  exhibit  a  moderately  high  Q  in  air 
(above  25)  -  For  these  (and  most  resonances  with  a  Q  above 
5)  the  program  with  set  frequency  ranges  for  data  taking 
works  well.  However,  these  same  transducers  may  have  very 
low  quality  factors  under  load  (in  water).  To  avoid  data 
review  over  a  frequency  range  beyond  the  upper  limits  of  the 
instruments  and  perhaps  overlapping  other  resonances  of  the 
transducer,  it  is  recommended  that  the  program  allowing  for 
operator  input  of  the  frequency  range  be  utilized. 

C.   CONCLUSIONS 

It  was  noted  that  most  transducers  have  more  than  one 
resonance.  To  be  a  useful  analytical  tool,  the  operator  of 
this  system  needs  to  know  the  primary  design  frequency  of 
operation  for  the  test  transducer.  Should  the  resonance  at 
this  frequency  cause  the  admittance  or  impedance  to  exceed 
the  limits  cf  the  Dranetz  meter,  the  operator  may  opt  for 
data  taking  of  the  type  to  fall  within  limits. 

Manual  collection  of  impedance  and  admittance  data  for  a 
single  transducer  in  two  aeiia  (air  and  water)  and 
calculations  similar  to  those  done  by  the  computerized 
system  utilizing  data  at  thrse  hundred  frequencies  requires 
on  the   order  of   fifteen  hours.     This  system   carries  ou: 
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these  operations  in  seventy-four  minutes.  For  most 
transducers,  only  on:  type  of  dita  (admittance  or  impedance) 
is  normally  collect ea.  A  typical  measurement  which  produces 
results  and  the  plots  similar  to  Figures  3.15  through  3.23 
would  take  fifty-two  minutes.  (Since  each  collection  of 
data  allows  for  four  plots  and  two  lists  of  data,  the  time 
to  run  the  program  collecting  impedance  and  admittance  data 
in  both  media  could  be  redaoed  by  forty  minutes  and 
collection  of  one  type  of  lata  in  both  media  by  twenty 
minutes  if  no  plots  or  lists  ware  made  ani  only  the  listing 
of  transducer  parametars  similar  to  Figure  3.23  desired.) 

The  major  time  saving  factor:  for  this  system  would  be 
found  through  the  replacement  of  the  Dranetz  Impedance  Meter 
with  an  instrument  that  does  aot  require  as  much  operator 
interface.  However,  semi-automation  as  provided  in  this 
system  constitutes  a  major  tine-saving  over  manual  means 
even  when  the  equipment  set-up  is  not  altered  in  a  major 
way. 

The  plots  produced  using  an  X-Y  plotter  and  sweep 
oscillator  enable  meisurements  of  frequency  with  an  accuracy 
of  about  one  percent.  Using  tha  computer,  the  precision  is 
better   than   one   ia    ten   tiousand.      Similarly,    for 
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calculations  using  real  and  imaginary  components  of  the 
impedance  or  admittance,  computer  precisicn  is  better  than 
one  in  ten  thousand  as  opposed  to  about  one  in  ten  for 
manual  manipulations.  However,  the  measurement  accuracy  is 
at  best  plus  or  minis  two  percent  for  the  Dranetz  meter. 
The  advantage  gained  over  manual  methods  by  using  the 
computerized  system  is  due  to  the  systematic  pauses  for 
Dranetz  meter  stablization  prion  to  all  measurements. 

This  system  has  been  designed  to  analyze  any  type  of 
transducer.  Specific  adaptation  to  hanile  a  particular 
transducer  type  with  known  resonance  range  can  drastically 
reduce  much  of  tfca  operator  interface  necessary  for 
evaluation  of  the  properties  of  the  transducer  and  reduce 
the  time  requirement  iccoriingly.  Adaptation  for  a 
particular  type  of  transducer  will  allow  variables  currently 
in  the  computer  program  to  be  made  constant.  Careful 
selection  cf  these  values  will  maximize  accuracy  and 
precision  in  the  desired  measurenents  and  calculations. 

D.   RECOMMENDATIONS  FDR  FURTHER  IMPROVEMENTS 

In  all  of  the  preceding  discussions,  it  must  be  apparent 
that  the  "weak  link"  in  this  system  is  the  Dranetz  impeiance 
meter  which   is  over  fifteen   years  old.     The  difficulties 
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wirh  this  device  ace  primarily  its  lack  of  computer 
controllability,  its  plus  or  minis  two  percent  accuracy,  and 
its  high  noise  at  low  signal  levels. 

Although  it  was  unfsrtuaat ely  not  available  for  this 
study,  the  Hewlett-Packard  4192A  Low  Frequency  Impedance 
Analyzer  seems  to  be  the  proper  replacement  for  the  Dranetz. 
It  is  capable  of  fill  program  control  via  the  HP-IB  and 
would  also  eliminate  the  nee!  for  a  ssparate  frequency 
synthesizer  and  data  acquisition  system.  It  is  highly 
recommended  that  this  device  bs  substituted  in  any  future 
applications  of  this  technique. 
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APPENDIX  A 
Two  programs  wers  written  to  implement  the  measurement 
procedures  and  the  calculations  to  evaluate  any  type  of 
transducer.  The  first  program  titled  "HYDRA2"  automatically 
looks  at  a  freguency  range  of  ten  times  the  bandwidth  on 
either  side  of  the  rssonance  frequency.  (This  is  modified 
for  water  to  be  five  times  ths  bandwidth  on  either  side.) 
This  program  also  automatically  selects  the  frequency  at 
which  to  take  shunt  measurements  based  on  the  Q  factors. 
The  second  program,  titled  "PICSBW",  allows  the  operator  to 
specify  the  bandwidth  for  anaLysis  and  the  frequency  at 
which  to  measure  shunt  values.  The  following  reference 
guide  is  applicable  for  both  programs.  Major  differences  in 
the  programs  occur  between  linss  2501  and  2650  and  between 
lines  8601  and  8635. 
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REFERENCE  GUIDE 

Line  Nos.  Description 

1-60  Declarations,  inputs,  "bookkeeping" 

60  -  90  Admittance  Spectrum  (freq.  &   Q)  using 

subroutine  7300 

100-149  Iapedance  Soectrum  (freq.  5  Q)  using 

subroutine  7300 

150  Obtains  blocked  values  using  subroutine 

2501 

153-176  Impedance   data    collection    (option) 

177-186  Admittance    data    collection    (option) 

198-450  Listing   and    plots   of   data    using 

subroutines 

451-500  Motional   data    calc.    for    magneto,     (for 

»*•    data) 

525-540  Motional   data   calc.    for    elec.     (for 

'P    data) 

541-642  List   of    motional    data   and    plot 

645-657  Motional   data    calc.    for    elec.     (for   'Z' 

data) 

658-668  Motional    data    calc.    for    magneto,     (for 

•Z1    data) 

670-757  List    of   motional    data   and    plot 

800-940  Calculations    for    elec.    coupling  in   air 

1000-1031  Calculations    for    magn.    coupling  in    air 

1040-1146  Calculations    for   elec.    coupling  in   water 

1150-1225  Zalculations    for    magn.    coupling  in    water 

1350-1410  Options   for    rerun   or   change   of   medium 

2501-2650  Subroutine    to    determine    blocked  values 

7000-7545  Subroutine    for    search   and    spectrum 

analysis 

8100-8110  Subroutine   to    print   data    list 

8200-8250  Subroutine    to    plot    real    da^a    vs.    freq. 

8300-8380  Subroutine    to    plot    imaginary    data    vs. 

freq. 

8400-8450  Subroutine    to    plot    'circles' 

8549-8573  Subroutine    to    obtain   normalization 

fa  ctor s 
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860  0-86  99 
8700-8799 
9000-9350 


Subroutine  to  collect  lata 

Subroatins  t3  determine  max/min  values 

Subroutinss  for  printouts 
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1  REM    **HYDRA2** 

2  OPTION    BASE     1 

3  SHORT   Q0,UlrQ8,Q9 ,w  ,E  ,E2  ,E3 , E5 ,R8 ,R7 ,X2 ,X3 ,X5 ,X6 ,X7 ,X8 
,G7,G8  f£2  ,£3,B5,B6  ,B7  ,B8  ,K,Kl,P5  ,P6  ,W 

7  SHORT   A,F2,F3,F4rF5,H3,H4,Jl,J2,J3#J4  ,R,X,G,B 

8  INTEGER   MfY#Yl,I,Y2fJ,Z2  ,H1,H2  £2, Si 

10  DIM    B$[2]  ,A$[2]  fTS.[lJ  ,C  (.200,3)  ,A1$.[4]  ,D$[30C]  ,L$[20]  , 
R$.[l]  ,I$[1]  ,S$.[1G]  ,C$[5J  ,C$.[2  3]  ,H$[2  3] 

11  PEAL    S  ,N  ,N8  ,N9  ,R0  ,R9,X0  ,X9  ,L  ,C0  ,L0  ,0  ,P  1  ,H8  ,H9  ,P,W0,G0 
,G9  ,B0  ,BS  ,N(  10)  ,M(10)  rF0  ,F1,F6  ,F8  ,F  9 

14  REAL    S6  ,E4,E6  ,L5  ,D  ,D1,D2  ,D3,D4 

15  Xb  ,D1,D2  ,D3  ,£4  =  0 

16  k,E2 ,E3,£5,R6 ,R7 ,X2,X3,X5=0 

17  Z3,A,F2,F3#F4#F5,B6,B7,B8=C 

18  S6.E4  ,E6,M,X,Y1 ,1  ,Y2,J  ,22=0 

19  KfKl,P5,P6,W,QO,Ql,g8,Q9fE=0 

20  CO,LO,0,P1,B,E8  ,H9  ,P  ,W0  ,S=0 

21  N  ,N6  ,N9  ,R0  ,R9  ,X0  ,X9  ,GG  ,G9=0 

22  £G,BS,H1,H2  ,F0  ,F1,F6  ,FS  ,F9=0 

25  DISP    ".THIS    PROGRAM   AUTOMATICALLY    SELECTS    HALFWIDTH    FO 
R    DATA    COLLECTION. " 

26  DISP    "THIS    LOOKS    AT    10    X    BW    IN    AIR    &    5    X    £W    IN    WATER. 

27  DISP    "TO    INPUT    YOUR    OWN    BANDWIDTH    USE  "PICKEW"    PFXGRA 
M    INSTEAD.  (HIT  "CCNT')" 

29  PAUSE 

30  CLEAR    §    BEEF 

.35  DISP  "ENTER  TYPE  OF  TRANSDUCER  TO  BE  MEASURED?  ("m" '  FO 
k  MAGNETIC  COUPLING  OR  "E  '  FOR  ELEC.)" 

42  INPUT  T$@  DISP  "TO  GET  A  COMPLETE  SET  OF  DATA  YOU  NEE 
D  MEASUREMENTS  IN  ECTH  AIR  AND  WATER.  " 

45  DISP  "DC  AIR  FIRST. IN  WHAT  MEDIUM  ARE  YOU  OPERATING?  ( 
ENTER  "!'  FOR  AIR  OR  "2"  FOR  WATER)" 

46  INPUT  M€  CLE/^R  <s    BEEP 
4  7  IF  M=2  THEN  4  9 

4  8  C$.  =  "AIR"  @  Z2  =  I0  %    GOTO  51 
49  CS=MWATER"  @  Z2  =  5 

51  PRINTER  IS  2  &  PLOTTER  IS  1 

55  DISP  "SET  THE  LRANETZ  FOR  ADMITTANCE  .S El  CN  LOWEST  SC 

ALE  THAT  wCN '.1  PEAK  DURING  FUN.iCCNl)" 

6C  PAUSE 

65  PRINT  £  PRINT  %    PRINT 

70  PRINT  "ADMITTANCE  IN  ";C$ 

75  GCSUE  7000 

8  0  IF  M*2  THEN  9  0 

05  F0=F6  @  w0=^  @  GOTO  100 
90  F1=F6  @  Q1*Q 

100  GCLEAR  @  CLEAR  £  BEEP 

105  DISP  "SET  DRANETZ  FOR  IMPEDANCE- (SCALE  NEEDED) -HIT  C 

ONT  " 
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110  PAUSE 

115  PRINT  @  PRINT  §  PRINT 

120  PRINT  "IMPEDANCE  IN  " ;C$ 

125  GOSOB  7000 

130  IF  tt*2  THEN  14  0 

135  F8=F6  @  C3=Q  @  GOTO  141 

140  F9=F6    @   C;9=<j 

141  CLEAR    <s    EEEP 

142  PRINTER    IS    701,76    £    PRINT    USING    144 

14  4    IMAGE    2/ 

145  PRINT  "F(Y-AIR)=  "  ;  INT  (F  C)  ,  ,"F  (  Y.-WA1)  =  "  ;  INI  (F  i )  ,  ,"  F 
(Z.-AIR)  =  ";INT(F8)  ,  ,"  F  (Z-WAT)  ="  ;  INT  (F  S) 

146  PRINT  "Q(Y-AIR)*  " ;Q0f  ,aQ ( Y-WATEE) *  " ;QL f ,MQ<Z-AIR) - 
";Q8  ,  ,"C(Z-wATER)  =  ";<*9 

147  PRINTER  IS  2 

148  CLEAR  £  DISP  "THESE  ARE  THE  INITIAL  ESTIMATE £.  (HIT  "' 
CCNT*'  WHEN  READ*  TO  PROCEED)" 

149  PAUSE 

150  GCSUE  250  1 
152  CLEAR  @  EEEP 

15  3  DISP  "TO  COLLECT  ADMITTANCE  DAI A  (FOR  ELEC.  COUPLING 
)  ENTER  "1'.  ENTER  "2"  FOR  IMPEDANCE" 

154  INPUT  L2(9  PRINTER  IS  2 

15  5  IF  L2=2  THEN  158 

156  RS  =  "G"  @  I$="E"  g  S$="ADMITTANCE"  d    H$="CCNDUCTANCE 
(MICRGMHOS)"  £  G$="S0SCEPTANCE  (MICRCMHGS)" 

157  GOTO  159 

158  RS="R"     £    I$="X"     d    S$=" IMPEDANCE"     £    G$=" REACTANCE     (CH 
MS)"     6    H$="RESISTANCE     (OHMS)" 

159  IF    L2=l    THEN    176 

160  CLEAR    g    DISP    "SET    ON    Z    AND    ENTER    SCALE    FACTCF    TO    COL 
LECT    DATA  .  " 

161  INPUT    S6@    PRINT    S  $ ;  "    IN    ";C$    C°    PRINT    "SCALE    =";S3;" 
OHMS" 

162  CLEAR    £    DISP    "I    AM    IvGRKING" 
164    GCSUE    8600 

16  5    IF    Y=2    THEN    152 

166  IF    M=2    THEN    16 S 

167  R8  =  E    (?    X3=E2    t    X6=E3    t    F4=E4 

16c    X6=£5    $    F5  =  E6    <s    F6  =  E1    &    GOTO    19b 

169  R7  =  E    £    X2  =  E2    (y    X7  =  E3    §    F4=£4 

170  X5=E5    i§    F5=£6    £    5  9=£l    g    GOTO     19  6 

176  CLEAR    &    BEEP 

177  DISP    "SET    CN    Y    S.    ENTER    SCALE    FACTOR    IN    MCFGMECS    TO 
COLLECT    LATA .  " 

176    INPUT    Sbt?    PRIN'I    S:?;"    It.    ";C£    &    PRINT    "SCALE    =";Sb;" 
MICRCMhOS" 

179  GOSLiE    6600 

180  IF    i=2    THEN    152 
lbl    IF    M=2    THEN    lb 5 
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162   G8=E    @    B3=E2    @    E8=E3    @    F2=E4 

363    B6=E5    S    F3=E6    (§    F0=E1    @   GCTC    198 

185  G7=E    g    B2=E2    G    E7=E3 

186  F2=E4    @    B5=E5    @    F3=E6    $    B 1=E 1 

198  PRINT    USING    199 

199  IMGE    2/ 

2CG  DISP  "DO  YOU  DESIRE  A  LIST  OF  COLLECTED  DATA  NEAR  RE 
SONANCE?  (100  POINTS)  (1  =  YES  l2=NC)" 

201  INPUT  Y 

202  IF  Y=2    THEN  300 
204  IF  L2=2  THEN  240 

210  PRINT  "ADMITTANCE  DATA  FOR  ";C$;"  IN  MICROMHCS"  tf  PR 
INT  USING  215 

215  IMAGE  2/ ," FREQUENCY"  ,  4X  ,  "PEAL"  f6X , "IMAGINAFY" 

216  GOTO  26  5 

240  PRINT  "IMPEDANCE  DATA  FOR  ";C$;"  IN  OHMS"  @  PRINT  US 
ING  215 

265  GOSUB  8100 

300  CLEAR  (9  DISP  "DC  YOU  DES1FE    A  PLOT  OF  DATA?  (  1  =  Y  IS  ,  2  = 
NO )  " 

30  5  INPUT  Y@  GCLEAR  @  CLEAR 
315  IF  Y=2    THEN  450 

.320  DISP  "ENTER:  1=  G/R  VS.  FRE^;  2=  B/X  VS.  FREQ;  3=  B 
VS.  G/X  VS.  R;  4=  END  PLOTTING  LOOP." 
32  5  INPUT  Y 
330  IF  Y=2    THEN  385 
3  35  IF  X=3  THEN  4  25 
340  IF  Y  =  4  THEN  450 
34  5  GOSUB  8  20  0 
356  LDIR  0  @  PEN  1  (y  PENUP 

3  57  MOVE  D( 1,1) +.2*  (D (300,  1) -D  (  1,1)  )  ,-(E*.6) 
36 C  LABEL  R$;"  VS.  FREQUENCY  IN  ";C$ 
361  GOTO  320 
':8  5  GOSUB  b30  0 
3  97  LDIP  0  t  PEN  1  (§  PENUP 
398  MOVE  D(l,l)+  .2*  |D(.20G,1)-D  (  1  ,1)  )  ,-L-.5*L 

400  LAEEL  1$;"  VS.  FREQUENCY  IN  "  ;  C  $ 

401  GOTO  32  0 

425  GCLEAR  &  CLEAR  y  GCSUB  t4UG 
431  PEN  1 

442  LDIR  C  £  PENUP 

443  MOVE  -(.5*L)  r-L-.2*L 

445  LAtEL    "INPUT    ELECTRICAL    "  ;  S  § ;  "    PLOT     ( "  ;I  $;  M  vS .  M  ;  £$?  " 
)     IN    ";C$ 

446  GOTO    320 

450  IF    L2=2    THEN    645 

451  IF    T$="E"    THEM    52  5 

h52    GCLEAR     fc    CLEAR    @    DISP    "I    AM    CALCULATING    MOTIONAL    CA1 

A" 

465    C0=E* .000001/ (2*PI*t6) 
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466  LO=l/( (2*PI*F6)  2*C0) 

470  FOR  1=1  TO  300 

475  D(I  ,2)=D(I  ,2).-G 

480  D(I  ,3)=D(I  ,3).-D(I  ,1)*B/F6 

48  5  NEXT  I 

450  GOSUB  8700 

500  GOTO  541 

525  GCLEAR  £  CLEAR  <s    DISP  "I  AW  CALCULATING  MOTIONAL  DAT 

A" 

532  C0=B*.000001/(2*PI*F) 

535  FOR  1=1  TO  300 

536  D  (I  ,2)=D(I  ,2).-G 

537  D(I,3)=D(I  ,3).-D(I,l)*S/F 

5  38  NEAT  I 

540  GOSUB  8700 

54  1  CLEAR  @  GCLEAR  £  BEEP 

542  DISP  "DC  YOG  DESIRE  A  LIST  OF  MOTIONAL  DATA  NEAR  RES 
ONANCE? (1=YES,  2=NC)" 

543  INPUT,  Y 

545  IF  v=2  THEN  604 

575  PRINT  "MOTIONAL  ADMTTANCE  DATA  IN  ";C£  (?  PRINT  USIN 

G  580 

58  0  IMAGE  2/,  "FREQUENCY"  ,  4X  /'REAL"  ,6X,  "IMAGINARY" 

595  GCLEAR  £  CLEAR 

600  GOSUB  8100 

604  CLEAR  £  SEEP  @  BEEP 

60  5  DISP  "WAN!  A  PLOT? (  1=YES  ,2=NC) " 

610  INPUT  Y2 

bl5  IF  i2=2  THEN  755 

bl6  PLOTTER  IS  70  5  §  PEN  1 

620  GOSUE  8400 

628  LDIR  0  £  PENUP  (§  PEN  1 

629  MOVE  - (. 5*L) ,-L-.2*L 

630  LAEEL  "MOTIONAL  ";S$;"  PLOT  FOR  ";C$ 
642  GOTO  755 

6  45  CLEAR  fc  GCLEAR  @  DISP  "I  AM  CALCULATING  MO'ilCNAL  DAI 
A" 

646  IF  T$="M"  THEN  65o 
64  7  FOR  1=1  TO  30  0 

648  D(I #2)=D( I ,2)-R 

649  C0=1/(2*PI*F*X) 

650  L0=l/( (2*PI*F6) ~2*CG) 

651  D  (I  ,3)=D(I  ,3).-2*PI*D(I  ,1)  *L0 
b52  NEXI  I 

653  GOSUB  8700 

654  C0  =  AES(CQ)  <?  L0=ABS(L0)  g  BEEP 

6  35  DISP  "DC  YOU  /.ANT  A  LIST  CF  THE  CAIA  NEAP  RESONANCE? 
(  1=YES  ,2  =  NO)  " 
6  56  INPUT  Y 
657  GOTO  670 
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658    L0=X/(2*PI*F6) 

662  FOP    1  =  1    TC  .200 

663  D(I ,2)=D(I  .2)  -R 

664  D(I  ,3)=D(I ,  3)-D(I , 1)*2*PI*LG 

665  NEXT    I 

666  GCSUB    8 7 CO 

66  7    BEEP    @    CISP    "    DC    YOU    VvANT   A    LIST    OF    CATA   NEAR    RESONA 
NCE?(1=YES ,2=NQ)" 
668    INPUT    Y 

6  70    IF    *=2    THEN    706 

66  5  PRINT  "MOTIONAL  IMPEDANCE  CATA  IN  ";C$  @  PRINT  USING 

666 
686  IMAGE  2/, " FREQUENCY ", 4X  ," REAL"  ,6X ," IMAGINARY" 

7  00  GCLEAR  &  CLEAR 

705  GCSUB  8100 

706  CLEAR  @  BEEP  @  BEEP 

710  DISP  "WANT  A  PLOT? ( 1=YES , 2=NO ) " 
715  INPUT  Y2 

720  IF  Y2=2  THEN  755 

721  PLOTTER  IS  70  5  @  PEN  1 
725  GCSUB  6400 

737  LCIR  0  e  PENUP 

736  MOVE  -  (. 5*L)  ,-L-.2*L 

740  LABEL  "MOTIONAL  ";S$;"  PLOT  IN  ";C$ 

75  5  GCLEAR  @  CLEAR 

756  PRIN1  USING  757 

757  IMAGE  3/ 

800  DISP  "I  AM  DOING  CALCULATIONS  FOE  XOU.  PLEASE  o£  PAT 

IENT." 

60  3  IF  M=2  THEN  104G 

810    IF    T$="M"    THEN    1000 

8  25    IF    L2=2    THEN    6  60 
640   Q0*E1/AES (E4-E6) 

o45    Cl=A£S  (E2/ (2*PI*FG.)  ) 

6  55    GOTO    68  5 

860    w8=El/ABS (E4-E6) 

o65    C1=A£S(1/ (2*PI*FS*E2) ) 

685    K=  i-(F0/Fb) "2 

695  Kl=Cl/(AbS (CC)+C1) 

5  30  GCSUB  9  000 

9  21  Pf INTER  IS  2 
94  0  GCTC  12  50 
1000  K=1-(F6/F0) "2 
lu05  IF  L2=2  THEN  lo25 
1010  L,G=El/AES(E<i-E6) 
L015  GCTC  10  30 

1025  ^b=El/A£S (E4-E6) 

1030  GCSUB  920o 

1031  PRIN'iER  IS  2  e  GOTO  1250 
lu40  IF  T$*"M"  THEN  1150 

98 


104 1  DISP  "If  YGU  HAVE  DATA  FOR  AIR,  ENTER  1;  ELSE  2" 

1U4  2  INPUT  Y 

104  3  RAD 

1U44  IF  L2=l  THEN  1110 

1045  IF  Y=2  THEN  1062 

1060  P6=D4*  (D3-D4)/(D3*R7) 

1061  W0»F1 

1062  Q1=E1/AB£(E4-E6)  (2  K*1-(F  1/F9)  *2 
1070  GOSUB  9100 

1076  PRINTER  IS  2  ^  GOTO  1350 
1110  IF  Y=2  THEN  1135 

1126  wc=fi 

1130  P6=D2*  (D1-D2)/(D1*G7) 
1135  gl=El/AES(E4-E6) 
1140  K  =  1-(F1/F9)  "2 

1145  GCSUB  9100 

1146  PRINTER  IS  2  §  GGTC  1350 

1150  RAD 

1151  DISP  "IF  YGU  HAVE  DATA  FCR  AIR,  EN1ER  1;  ELSE  2" 

1152  INPUT  Y 

1153  IF  L2=l  THEN  1160 

1154  IF  Y=2  THEN  1162 

1155  B=.5*ACS(  (R7-R9)/D4) 

1156  Hb=SIN(E)~2  V   H9=COS(E)~2 

115  7  P5=(  (1+D3/K0*H9)  *.  5-  (L-D3/R0*H8  )  ~.  5  )/  (  (1+D3/R0*H9) 
.5+  (1-D3/R0*H8)  ".5) 

1159  F=D4*SIN (2*B)/(4*R9*£S)  @  vvG=F9*  (  2*P  +  (  P  ~2+  1)  *, .  5) 

1160  P6=D4* (D3-D4)/(D3*R7) 

1162  C9=E1/ABS(£4-E6)  @  K= 1-(F 9/F 1) "2 

1163  PRIN.'I  "£=";£, "K=";K,"Kl=";Kl,"H£  =  "  ;ti8  ,"H9="  ;H9 

1164  GGSUE  9300 

1168  PRINTER  IS  2  &  GGTO  1350 
1180  RAD 

116  2  IF  Y=2  THEN  1210 
1184  B=.5*ACS(  (G7-G9.)/D2) 

1186  H8  =  SIN(S)~2  @  h9=CCS (B) "2 

li£8    P5=(  (i+Dl/GG*H9)  *.  5-(  1-D1/GG*H8)  \  5)/  (  ( 1+Dl/GU*H&) 

.5+  (1-Dl/GG*H8) * .5) 

1190    P«*D2*SIN  (2*B)/(4*G9*C.U     i    w0*Fl*  (2*P+(P*2+1)  "  .5) 

L195    Pb=C2* (D1-D2)/  (Dl*G7) 

1210  QS*E1/ABS(E4-E6)     @    K*1-(F9/F1)  *2 

1211  PRINT    "E  =  ";E,"K=";Kr"rvl  =  "  ;K1,"H6  =  "  ;Ho  ,"n9=";H9 
1215    GGS06    9300 

1225    PRINTER    IS    2    €    CLEAR    Eg    EEEP 

135U    DISP    "DC    YGU   DESIRE   AKCTHER    IN    THIS    MEEIUM?(TG    GET 

THE    OTHER    TYPE    DATA?) ( 1*YES ,2*NO) M 

13  55    INPul    Zl 

13  57  PLOTTER  IS  1 

1360  IF  Zl*l  THEN  152 

1361  PRINT  USING  126  2 
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1262  IMAGE  3/ 

1365  DISP  "DC  YOU  DESIRE  A  RUN  IN  ANOTHER  MEDIUM? ( 1= YES , 

2=NG)  " 

1370  INPUT  Zl 

1375  IF  Zl=l  THEN  45 

1360  DISP  "JUST  FOR  THE  RECCRC  ,  INPUT  WATER  TEMP,  AIR  TE 

MP,  TRANSDUCER  SER  NG ,  MODEL  NO.," 

1381  DISP  "AND  TYPE  (E  OR  M)" 

1385  DISP  "INPUT  "U"  IF  INFO  IS  UNKNOWN" 

1390  INPUT  W  ,A  ,S  ,N  ,L$ 

1391  PRINTER  IS  701  ,76 

1400  PRINT  "wATER  TEMP=" ; w  ,  ,"AI R  TEMP*" ;A  ,  ," SER.  NG.  =  ";S 

,  ,"  MODEL  NO  .  =  "  ;N  ,  ,"  TYPE  "  ;  L$ 

1410  GOTO  S998 

2501  REM  *SHUNT* 

250  3  CLEAR  %   BEEP  @  DISP  "WE  ARE  FINDING  SHUNT  VALUES  FC 

R  G/£  OR  R/X." 

2510  A$="FR"  £  E»="HZ"  &  Sl=l 

2511  IF  T$*"M"  THEN  2520 

2512  IF  M=2  THEN  2516 

2513  Q=MAX(^0,*6)  @  F=F0- (Z 2-2) *F0/w 

2514  IF  F>2C0  THEN  2535 

2515  F=2GG+(FG-2*FG/C.-2CC)/5G  @  GOTO  2535 

2516  Cj=MAX  (Cl,w9.)  (s  F=F  3-  (Z2-2) *F1/Q 

2517  IF  F>200  THEN  2535 

251c  F=20U+(Fl-2*Fl/w-20G)/5G  @  GOTO  2535 
25  20  IF  M«2  THEN  2  52  4 

2521  Q«MAX(W0#Q8)  £  F=F6-  (Z  2-2)  *F6/^ 

2522  IF  P>20G  THEN  2535 

2523  f=200+(Fa-2*Fb/C-20G)/50  @  GOTO  2535 

2524  fi*MAX(Sl,gS.)  g  f  =t  9-  (Z  2-2)  *fS/w 
25  25  IF  F>20  0  THEN  2  5  35 

2527  F  =  2CC+(f9-2*r  9/t*-2CG)/ 50 

2535  D$=vAL$(F) 

25  3  7    OUTPUT    717    ?A$rD$,B$ 

2  5  36    WAIT    10 CO 

25  40  DISP  "SET  DRANETZ  FREw .  SCALE  10  CCvER  FREg .  CN  SYN 

THESIZER  .ZERO  METERS  .SET  ON  FS.(CONT)" 

2541  PAuSE 

25^2  DISP  "HERE  aE  GE1  THE  NORMALIZATION  FACTORS." 

2543  OUTPUT  70^  ;"  vll" 

2546  GOSUB  b549 

2550  G  ,B,R,X=0 

2551  A^  =  "FR"     g    E$.  =  "HZ"     @    D$=VAL$(F) 
2553    CUTPUI    717    ;A$,D$,B;? 

2561  CLEAR    @    DISP    "SET    DFANETZ    ON    i  ,SET    SCALE  ,S El    CN    NCF 
M   EOF    FILTER    +    PHASE,PLUG    IN    TRANSDUCER." 

2562  aEBP    i    DISP    "ENTER    SCALE    FACTOR    IN    MICRCMHOS.SAME    5 
CALE   AS    FOR    DATA.(hE    GET    SHUNT    G/£)  " 

25c4    INPU1    S6@   CLEAR    fc    DISP    "I    AM    WORKING    rO    GE1    SHUNT    i/ 
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ALUES"  (s  BEEP 

2565  FOR  1=1  TO  10 

2566  OUTPUT  717  ?A$,D$,£$ 

2567  WAIT  1C0G 

2566  OUTPUT  709  ;"AC3" 
2  569  ENTER  709  ;  M(I) 
257  0  G  =  G+iy.  (1) 
2571  OUTPUT  7U9  ;  "AC4" 


2  57  2  ENTER  7  09  ; 
257  3  B=£+N(L) 

2574  NEXT  I 

2575  G=G/(I-1)  i 

2576  G=G*S8/N6  @ 


N(I) 


E  =  £/(I-]) 
B*B*S8/N9  fc  BEEP 

2578  DISP  "SET  LRANETZ  ON  Z ;  SET  SCALE  FOR  NAX  RESPONSE. i^ 
CVE  TRANS.  INPUT.  ENTER  SCALE  fACTCR" 

2579  INPUT  Sc(y  CLEAR  @  DISF  "I  A*  wOFKING  TG  GET  SHUNT  V 
ALUES" 

2560  FOR  1=1  1C  10 

2581  OUTPUT  717  ?A$,D$,B$ 

2582  WAIT  1000 

258  3  OUTPUT  709  ;"AC3" 
2  564  ENTEF  709  ;  b:(l) 
2565  R=R+M(I) 


;"AC4" 

N(I) 


25c6  C0TPLT  709 
2  58  7  ENTER  709 
2566  X=N(I)+X 
25b 9  NEXT  I 

2590  R=F/(I-1)  &  X=A/(I-]) 

8/N  9 


2591  F=R*S8/N8  fc  X 


=  x*.c 


<.592 
25^5 
25. so 
2  6oL 
2601 
2  60  2 
2  60  3 
2  60  5 


IF  T$="M°  THEN  2b  10 
IF  N=l  THEN  26 Ol 
R9=R  <e  X9=X  fc  G9=G 
G0=C  |  E0=E  &  R0=R  & 
PRINT  "SHUNT  DALLES" 
PRINT  USING  26  0  3 
INAGE  1/ 
PRINT  "G0=", -GO  ,  ,"E0  =  " 


BS*B 

Xo  =  X 


It    GOTO    26  01 


;GS,  ,"E9=";E9  ,  ,"  R9="  ;R  9  ,  ,"  X9=  "  ;  X9 
26C6    PR INI    6   GCTO    26  3 u 
261C    IF    Sl=2    THEN    2620 

2611  J  1=G    £    0  2=8    (9    0  3  =  R    t    J  4  =  X 

2612  IF  N  =  2  THEN  2ol5 

2613  F=F6+(Z2-2) *F6/g 
2615  F  =  F9+ (Z2-2)  *F9/i* 
2620  G=.5*fG+01)  £  ii=.5*(S+J2) 
X)  i    GC'IC  2595 
2630  GCLEAR  %    CLEAR  @  6EEF 
2c  3  5  EISP  "EN1ER  "l"  IF  ALL 
EPEAT." 
2640  INFOT  Y 


L  ,   ,"R0=";R0  ,  /' X0  =  "  ;X0  ,  ,"G9  = 


Sl*2 

Sl=2 


GOTO 

GGIC 
t  R  = 


2550 
2550 
5*  (R+J3) 


X*.5* (J4+ 


IS  sitLL 


IF  /Co  NEED  A  B 
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2  6  45    IF    Y=2    THEN    2  501 

2  650    RETURN 

70CG    REM    **F-Q** 

7001    GCLEAR    @    CLEAR 

70C5   DISP    "ENTER    LOWER  ,    UPPER    FREQ.    FOR    SWEEP. IF    RESCNAN 

CE    UNKNOWN    USE    10000,200000.     (F1,F2)" 

7010  INPUT    F4,F2 

7011  DISP  "ENTER  BhS    VOLTAGE  NEEDED.  (3000  MV  FOR  LOW  FRE 
Q.  AS  LOW  AS  2000  MV  FOR  HIGH)" 

7012  INPUT  A£  Al$=VAL$(A) 
7o20  CLEAR  g  BEEP 

7025  A$="AK"  @  B$="NR" 

7035  OLTFUT  717  ;A$,A1$,B$ 

7040  OUTPUT  709  ;"AC3vT3" 

7045  F3=CEIL( (F2-F4)/3GC) 

70b  C  CLEAR  @  DISP  "I  AN  WORKING  TO  GET  THE  SPECTRUM  DATA 

FOR  YOU,lHEN  WILL  SHOW  &  MAKE  A  PLOT" 
7065  FCR  1=1  TO  300 
7070  D  (I  ,1)=F4+I*F3 
7075  A$="HZ"  £  B$.  =  "FR" 
7063  D$=VAL$(D( I  ,1)  ) 
7  08  5  OUTPUT  717  ;B$,D$,A$ 
7C9C  OUTPUT  709  ;"AC3VT3" 
7095  wAIT  100 
7100  ENTER  70  9  ;  D( I  ,2) 
710  5  NEXT  I 

7110  F5  =  F4-.1*(F2-F4)  (?  F6=F  2+  .  1*  (F  2-r  4) 
7115  F6=F2-h.l*  (F2-F4) 
7 US  PLOTTER  IS  1  t  GCLEAR  £  CLEAR 
7125  SCALE  F5  ,F 6  ,- . 1  ,1 . 2 
7130  XAXIS  0  ,2000  ,F4  ,F2 
7135  XAXIS  F4f  .1  ,-.1  ,1.2 
714  0  FOR  1=1  TO  30 C 

714  5  PfiNUP 

7150  PLOT  D(I  ,1)  ,D( I  ,2) 

715  5  NEXT  I 

7156  GRAPH  @  CCFY 

7160  DISP  "NEED  ANOTHER  FUN?  (  1= YE S ,2=NC) " 

7161  INPUT  1 

lit!    IF  i*l  THEN  7  0u5 

716  3  CLEAR  %    GCLEAR  g  BEEP 

71b5  DISP  "ENTER  DECISION  POINT  FOR  AMPLITUDE  (C  TO  1.2) 

.  {LESS    THAN  THE  MAX  DISPLAYED  ." 

7170  INPUT  S@  CLEAR  @  BEEP 

7160  PRINT  "AMPLITUDE  IN  VOLTS";"    FREQUENCY" 

7181  PRINT 

719  0  FOR  1  =  1  TO  30  C 

llbi    IF  D  (I  ,2)  <S  THEN  7210 

7  2UU  PRINT  USING  7  20  5  ;  D ( I ,2 ) ,D ( I , 1 ) 

7205  IMAGE  IX  ,C . LDCDCEDL , 1  OX  .DDODDC  .DD 
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7  2 1C  NEXT  I 

7  215  PRINT  USING  7  22  0 

7220  IMAGE  3/ 

7225  CLEAR  §  BEEP 

7226  DISP  "IF  YOU  ARE  FEADY  TC  CONTINUE  ENTEB  "l";  OTHER 
WISE  ENTER  '2"  ." 

7227  INPUT  Y 

7228  IF  Y  =  2  THEN  716C 

72  35  DISP  "WHAT  IS  CENTER  FREQ,  HALFWID1H  TO  CONSIDER?" 
7236  INPUT  F,N(s  CLEAR  (s  BEEP 

7240  DISP  "I  AM  WORKING  TO  GET  GOOD  FREQ    AND  Q  DATA" 

7241  FOR  1=1  TO  30u 

7242  D(I  ,  1)=F.-N+I*N/15C 

7243  A$="HZ"  @  B$="FR" 
"/245  D$=VAL$(D  ll  ,1)  ) 
7255  OUTPUT  717  ;B$,D$,A$ 
7260  OUTPUT  70S  ;"AC3VT3" 
72b5  WAIT  100 

7270  EN1ER  70S  ;  D(I,2) 
7  275  NEXT  I 

7  28  5  CLEAR  @  GCLEAR  @  DISP  "I  AM  FINDING  THE  ABSOLUTE  MA 
i  AND  FREQ  UPPER  AND  LOWER" 

290  B1,B4=0  @  H,Hl  ,H2  ,H3,H4=1 

2S5  E2,B3=50 

300  FOR  1=2  TC  300 

310  IF  D(I  ,2)  <D  (H  ,2)  THEN  "/340 

315  IF  D(I  ,2) =D (H ,2)    THEN  7  325 

325  A6=D(I,2)  £  F6  =  D(I,1)  @  H=I 

3  40  NEXT  I 

345  A7=A6/2 

3  50  FOR  1=1  TC  H 

355  IF  A7=D(I,2)  THEN  7405 

360  IF  A7<D(I,2)  THEN  7365 

365  IF  D(I,2)<E1  THEN  7425 

370  B1=D(I,2)  g  H1  =  I  y  GOTO  7<,2  5 

36  5  IF  D(I,2)>B2  THEN  7400 

390  B2=D(I ,2)     @  H2=I 

4CU  GOTO  7  4  25 

405  Bl,B2  =  D<I,2)  9  Hl,H2=I  $  F7=D(H2,])  @  GOTO  7440 
7  42  5  NEXT  I 

7430  X=(A7-B1)/(B2-B1) 
74  3  5  F7=X* (D(H2,1)-D(H1,1))+E  (HI  ,1) 
744U  FCR  I=H  10  300 
74  4  5  IF  A7=D(I,2)  THEN  7  <i  S- 5 
7450  IF  A7>D(I,2)  THEN  7475 
7455  IF  D(I,2)>B3  THEN  7470 
7460  B3=D  (I  ,2)     <t    H3=I 
7^70  GOTO  7  515 

7475  IF  D(I,2)<=B4  THEN  7^50 
7460  B4=D ( I  ,2)     %    H4  =  I 
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7490  GOTO  7515 

7495  B3,34=D(I,2)  @  H3,H4«I 

7505  L5=D(H3,1) 

7510  GOTO  7530 

7515  NEXT  I 

7520  X=(A7-B4)/(E3-E4) 

7  525  L5=-(X* (D(H4  ,1).-D  (H3  , 1).)  )+D(H4  ,1) 

7530  C=F6/(L5-F7) 

7535  PRINT  "CENTER  PREQ  IS  ";F6 

75  36  PRINT  "Q    IS  ";C. 

7  536  DISP  "IF  YOU  ARE  READY   TO  PROCEED,  ENTER  ''!'.  TO  RE 

RUN  FOR  BETTER  VALUES,  ENTER  "2"'." 

75  3  9  INPUT  Y 

7540  IF  Y=2  THEN  7235 

7545  RETURN 

8100  REM  **DATA  LIST** 

8101  PRINTER  IS  2 

6102  FOR  1=99  TC  199 

6103  PRINT    USING    8105    ;    D(I  ,1)   ,D(I,2)  ,D(I,3) 

8105  IMAGE    DDDDDD.D,2X  .D.DCDE        ,2X,D.DDDE 

8106  NEXT    I 

8107  PRINT    USING    8106 
810b    IMAGE    3/ 

8110  RETURN 

8200  REM  **PLOT  RE** 

6201  GCLEAR  £  CLEAR 

8202  H8=D(1,1)-. 15* (D(300  ,1)-C (  1,1)  ) 

620  3  H9=D(300,l)+.l* (D ( 300 , 1) -C ( 1 , 1 ) ) 

62C5  PLOTTER  IS  705  (y  PEN  1 

8209  SCALE  H8  ,H9  ,- ( . 8*E )  , E  +  . 25*E 

6210    XAXIS    0  ,(D(  300  ,1)-C  (  1,1)  )/10  ,C  (  1  ,  1)  ,0  (.300  ,1) 

8  211    YAXIS    0(1,1)  ,E/5,-(.S*E)  ,E  +  .2*E 

6212    PENUP 

6215    FOR    1=1    TO    300 

8  216    PLOT    D(I,1)  ,D( I  ,2) 

8217   wEXT    I 

8216    PENUP    <s    DEG    $    LDIR    0,SIN(90) 

6221    FOR    L1=D(1,1)    TO    0(200,1)     STEP     (  D  (  30  0  ,  1 ) -O  ( 1 ,  1 )  )  /  lb 

!§    PENUP 
b222   MOVE    LI,- (. lb*£) 
6  22  3    LABEL    INT (Ll) 
6224    NEX'l    Ll 

6226  LDIR    0    @    PENUP 

6227  MOVE    D ( 1 , 1)  +  .2* (D ( 300  , 1) -D  (  1  ,  1 )  )  ,- ( . 7*E ) 

6228  LABEL  "FREQUENCY  (HZ)" 
8  2  30  PENUP  (9  PEN  1 

8232  LDIR  0 

b234    FOR    Ll=-(.b*E)     TC    E+.2*E    STEP    E/5    §    PEI.UP 

82  3  5   MOVE    D(  1,1) -.09*  (D(200,l) -0(1,1))  ,L1 

o2  36    LABEL    IN'l(Ll) 
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8237  NEXT  LI 

6240  LDIR  C,SIN(9C) 

8  241  MOVE  D(1,1)-.1*(D(3G0,1)-D(1,1)),.1*E 

824  3  LABEL  H$ 

6  2  50  RETURN 

8  300  RET-:  **PL01  im** 

6301  GCLEAR  (§  CLEAR 

8  30  3  H6=D(1I1).-.15*  <D(300,1)-D  (  1,1)  ) 

6  30  4  H9=D(3G0  ,1)+.1*(D  (300,1) -0(1,1)) 

8  305  L=MAX(AES(E5)  ,AtsS(E3)  ) 

6306  PLOTTER  IS  705  @  FEN  1 

8312  SCALE  H8  ,H9  ,-L- . 5*L  ,L+ . 25*L 

6313    XAXIS    0  ,(D(3G0,1)-C (1,1)  )/ lu  ,D(  1  ,  J)  ,D(  300  ,1) 

8314    YAXIS    C( 1,L)  ,L/5  ,-L-.2*L ,L+.2*L 

8  315    PENUP 

8  316   FOR    1=1    TC    30  0 

8325    PLOT    C(I  ,1)  ,D(I,3) 

8330  NEXT    I 

8331  PENUP    @    DEG    <e    LDIR    G,SIN(90) 

6333    FOR    L1  =  C(1,1)    TO    D(300,l)    STEP    (D  (  30  0  , 1 ) -D  (  1  ,  1)  )  /  10 

<d    PEN  OP 
8335  MOVE    LI  ,-L-.2*L 
6  3  36    LABEL    INT (LI) 
6337   NEXT   LI 
8340    LDIR    0    §    PENUP 

8  341   MOVE    D(1,1)  +  .2*(D(3GG,1)-D(1,1))  ,-L-.4*L 
8  34  2    LABEL    "FREQUENCY     (HZ)" 
6  345    PENUP    @    PEN     1    @    LLIB     0 

6352    FOR    Ll=-L-.2*L    TC    L+  .2*L    STEP    L/ 5    %    PENUP 
8355    MOVE    D ( 1 , 1) -. 09  *  (D (  300  , 1 ) -D  (  1  ,1 ) )  ,L1 
8  356    LABEL    INT (LI) 
8  3  57    NEXT    LI 

8360    PENUP    @    PEN    1    @    LDIR    C,SIN(S0) 
b 3  61   MOVE    D(1,1)-.1*(D(3GG,1)-0(1,1)),-(L*.5) 
6  366    LABEL    G$ 
o36  0    RETURN 
8  400    RE^.    **CIRCLES** 
64U1    GCLEAR    @   CLEAR 
8  4  02    L1=MAX(ABS  (E3)  ,AES  (E5)  ) 

6403  L=MAX(AES(E)  ,L1) 

6404  PLOTTER    IS    705    ts    PEN     1 

8406    SCALE    1 . 25* (-L- . 1*L)  ,  1 .  25*  (L+ . 3 *L )  ,-L-  .  2*L  ,L+ .  15*L 

6407    XAXIS    0  ,L/5  ,-L-.l*L  ,L+.  1*L 

6406    ¥AXIS    0  ,L/5  ,-L-.l*L  ,L+.  1*L 

6  409    PENOP    @    PEN     1 

84  11    FOR    1=1    TC    30  0 

6  412    PLOT    D(I,2)  ,L(I,3) 

64  13    NEXT    I 

64  14    PENOP 

b  4  15    LDIP    0    y    PEN    1    fs    PENOP 
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8  416  MOVE  - (L*.5) ,-L-. 1*L 

8  418  LAEEL  H$ 

6420  PEN UP  @  PEN  1  I  DEG 

8425  LDIR  0,SIN(9Q) 

8426  FCR  L1=-L-.1*L  TC  L+.1*L  STEP  L/ 5  @  PENUP 
6427  MOVE  Ll,-L 

8428  LAEEL  IN.T(Ll) 

8  42  9  NEXT  Ll 

6  435  LDIR  0  @  PEN  1  @  PENUP 

8436  FOR  Ll=-L-.1*L  TO  L+.1*L  STEP  L/5  §  PENUP 

6437  MOVE  - (.666*L) ,L1 

6436  LAEEL  INT(Ll) 

6439  NEXT  Ll 

8445  LDIR  0«SIN( 9G) 

8446  MOVE  -  (.95*L)  ,- (. 3*L) 
6446  LAEEL  G$ 

8450  RETURN 

8  5  49  REM  **NCRM** 

6550  N8  ,£m9=0 

b551  FOR  1  =  1  TC  1U 

6557  WAIT  5u0 

8  5  58  OUTPUT  7  09  ;"AC3" 

8559  ENTER  709  ;  N(I) 

8560  N6=N8+N(L) 

8565  OUTPUT  709  ;"AC4" 
6566  ENTER  709  ;  MI) 
8567  N9  =  N9+ML) 
b566  NEXT  I 
85b9  N8  =  N"6/  (1-1) 
8570  N9=-(NS/(I-1) ) 
8  57  2  PRINT  USINO  6  5  73 
6  57  3  IMAGE  3/ 

65  7  5  RETURN 

66  00  REM  ** VALUES** 
c6  0l  CLEAR  8  BEEP 

b6  02  DISP  "I  AMi  COLLECTING  THE  DATA  NOW" 

66  L  5  IF  L2=l  AND  M«2  THEN  66  2  0 

66  06    IF    L2=l    AND    M=l    THEN    6b 21 

66  07    IF    L2=2    AND    M=l    THEN    66  22 

66  6b    IF    L2=2    AND    M  =  2    THEN    o6  2  3 

6620    E6  =  F1    l§    si=lwl    (S    GCTC    bb25 

bbll    Fb=F0    @    ^=^0    @    GOTO    66  25 

8622    F6=F6    e   w  =  ^b    f    GOTO    b625 

66  2  3    F6  =  FS    g    <*=<*9 

c-625    FCR    1=1    1C    50 

o026    Lc    F6-Z2*F6/Q  +  (Z2-2)*F6/(50*vi)  >20G   THEN    &631 

6629    D(I  ,1)«200+(F6-2*F 6/^-2  00)  *I/50 

co 30    GOTO    a 6 32 

66  31   D  (I  ,1)»F6-2.2*F6/G.+I*  (22-2)  *F6/(50*CJ 

ob32    D$*VAL*  (DU,1)) 
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8633    OUTPUT    717    ;A$,D$,E$ 

86  34    WAIT    150 

66  3  5   OUTPUT    709    ;"AC3VT3" 

66  36    ENTER    709    ;    nti 

6637   D(I  ,2)=H6*S8/N8 

86  36  OUTPUT  709  ;"AC4VT3" 

8639  ENTER  709  ;  H9 

o640  D  (I  ,3)=H9*S6/N9 

6641  NEXT  I 

66  50  FOR  1=51  TO  2  50 

66  51  L(I  ,l)=F6-2*F6/^+(I-50)*2*F6/ (1G0*L) 

8652  D$=VAL$(D(I ,1) ) 

86  53  CUTPUT  717  ;A$,D$,B$ 

bo  54  /,AIT  250 

66  55   OUTPUT    709    ;"AC3VT3" 

66  56    ENTER    709    ;    H6 

8657    D(I  ,2) =H6*S6/N8 

86  58    OUTPUT    7  09    ;"AC4VT3" 

86  5  9    ENTER    7  09    ;    K9 

6660  D (I  ,3) =H9*S6/N9 

6661  NEXT    I 

8670    FOR    1=251    TO    300 

86  71    D(I ,1)=F6+2*F6/0+(I-2  50) * (Z 2-2) *F6/ ( 50*  C) 

6672   D$«VAL$(D(I  ,1)  ) 

o673    OUTPUT    717    ;Ai?,D$,E$ 

86  7  4    WAIT    150 

86  75  OUTPUT  709  ;"AC2VT3" 

66  76  ENTER  7U9  ;  H8 

b6  77  C(I  ,2)  =H8*S8/N6 

b6  78  OUTPUT  709  ;"AC4VT3" 

66  79  ENTER  709  ;  H9 

o660  D  (I  ,3) =H9*S6/N9 

66  61  NEaT  I 

6662  CLEAR    @    DISP    "I    AM    FINDING    MIN/KAX    VALUES    &    ASSOC. 
FREUS" 

86b3   GOSUB    o700 

6b84    IF    L2=2    THEfi    8687 

b6  6  5    PRINT    "GiV:AX  =  "  ;E  ,  ,"  £Gi*'X  =  "  ;  L  1  ,  ,"  £Gi-^  =  "  ;E2  ,  ,"u\A  =  ";l;  , 

t"FBiVX="  ;£4,  ,"tM="  ;E5,  r"FBMI=";£6 

6666    GOIG    6bb7 

66  6  7    PRINT    "RftAX*" ;E  ,  #"  FRi*;X=  "  ;  £  1  ,  ,"XF;«;X  =  "  ;E2  ,  ,"  aNX=";E3  , 

,"  FXMX=";E4,  ,"XNI=";£5  ,  ,"  FXML  =  "  ;  E6 

8697    DISP    "IF    ALL    IS    wELL  ,    ENTER"!'.     ENTER  ''2'    TC    RETAK 

E    DATA." 

o69b    INPUT    Y 

a 69 9    RET U FN 

6  7  0b    REK    **^AXNTN** 

b"/03    h  ,'n  1  ,H2=  lou 

6705    £=D(lub,2)     @    £1 ,£k  ,£6  =  L ( 1Gb ,1)     y    E2,E3  ,E5*D  ( lOu  ,  2) 

o7Uo    FOR    1=10  1    TC    2bo 
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8  710  IF  C  (I  ,2)  <D<H,2)  THEN  6  72  5 

6711  IF  C (I ,2)=C(H  ,2)  THEN  8720 

8720  E=D(I,2)  e  E1*D(I,1)  @  E2=£(I,3)  §  H  =  I 

6  7  25    IF    D  (1,3)  <D(Hl  ,3)    THEN    8740 

6  726    IF    C(I  ,3)  =D  (Hi  ,3)     THEN    67  35 

8735    E3=D(I,3)     @    E4=D(I,1)     g  H  1  =  1 

87*0    IF    D(I,3) >=C(H2,3)    THEN    6751 

6745    H2  =  I    @    E5  =  D(I#3)     @    E6=D(I,1) 

8751   NEXT    I 

6760    Jl=D(Hf2)     <e    J2=(D  (H+l,2)-D  (H-l  ,2)  )/2    @    J3=(D(H  +  1,2) 

+D(H-l,2).-2*ul)/2 

87b2    F6=-(J2/(2*J3) ) 

6  7o4    A=J  i  +  J2+Fb+J  3*F6~  2 

o770    F6=D  (H  ,1)+F6*  (D(H  ,1)  -C(H-1,  1)  ) 

8775    £1=F6    @    E=A 

6766    D=ABS(E3-E5) 

6  76  9    IF    M=2    THEN    87  S3 

6  790    IF    L2=l    THEN    8792 

6791    C3=D    @    GCTO    6796 

8  7  92    D1=E    @    GOTO    8  7  96 

6  79  3    IF    L2=l    THEN    8  79  5 

87S4    D4=C    g    GCTO    6796 

8795    D2=C 

o796    PRINT    "Dl="  ;  D  1  ,  ,"C2=";D2  ,  ,"D3  =  ";C3  ,  ,"D4  =  "  ;D4 

879  9    RETURN 

9  00  0    REM    **CALC-A** 

90  01    PRINTER    IS    7  0  1,76 

9002    PRIM    @    PRIN'l 

9o0  3    PRINT    "VALUES    MEASURED    IN    AIR    FROM    ";S^>;"    DATA" 

90 0 4    PRINI 

9011  PRINI  "ELOCKED  CAPACITANCE*  " ;A£S (Co )  ; "  FARAES  AI  " 
;F; "hZ" 

9012  PRINT  "ELECTRICAL  QUALITY  FACTOR  =  ";v0 

9013  PRINT  "MECHANICAL  QUALITY  FACTOR  =  ";^o 

9014  PRINT  "ELECTRICAL  RES^bA^CE    =";r0;"  HZ" 

9015  PRIN'l  "MECHANICAL  RESONANCE  =  ";:'&;"  HZ" 

90 Id  PRINI  "DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 

NT  =  "  ;K 

b017    PRINT    "STATIC    COUPLING   COEFFICIENT*    ";Kl 

9021    PRINT    "SLOCKED    RESISTANCE  =  ";  Ro;  "CHI- 3  " 

9031    PRINI    USING    9  0  32 

i?0  32    IMAGE    2/ 

9  0  5U    RETURN 

9  loO    REM    **CALC-w** 

9101  PRINTER    IS    701  ,76 

9102  PRINT    g    PRINI 

9  1o3  PRINT  "VALUES   FOR  MEASUREMENTS  IN  *AI£R  FRCM  "  ;  £  ^  ; 

"  DATA" 

9lo4  PRIN'l 

bill  PRINI  "ELCCKEl    CAPACITANCE  =  ";AES(CU); 
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";F;"HZ" 

9112  PRINT  "ELECTRICAL  QUALITY  FACTOR  =  " ; Q  ] 

9113  PRINT  "MECHANICAL  QUALITY  FACTOR  =  " ;  w  S 

9114  PRINT  "ELECTRICAL  RESONANCE  =  ";P1;"  HZ" 

9115  PRINT  "MECHANICAL  RESONANCE  =  ";F9;"  H2" 

9116  PRINT  "DYNAMIC  ELECTROMECHANICAL  CCUFLINC  COEFFICIE 
NT  =  ■  ;  K 

9120  PRINT  "MECHANICAL  POWER  UTILIZATION  FACTOR" ; AES ( P6  ) 

9121  PRINT  "  STATIC  COUPLING  COEFFICIENT*" ;K1 

9122  PRINT  "ELCCKED  RESISTANCE  =  "  ; R9; "OHMS " 
9132  PRINT  USING  913  3 

913  3  IMAGE  3/ 

9  150  RETURN 

9200  REM  **CALC-A/M** 

9201  PRINTER  IS  701  ,76 
92  02  PRINT  g  PRINT 

9203  PRINT  "VALUES  MEASURES  IN  AIR  FROM  ";S$;"  DATA" 
920  4  PRINT 

9211  PRINT  "BLOCKED  INDUCTANCE*?  ";ABS(LC);"  HENRIES  Ai  " 
;F6; "HZ" 

9212  PRINT  "ELECTRICAL  QUALITY  FACTOR  =  ";C.G 

9213  PRINT  "MECHANICAL  QUALITY  FACTOR  =  ";Q8 

9214  PPIN.T  "ELECTRICAL  RESONANCE  =";F0;"  HZ" 

9215  PRINT  "MECHANICAL  RESONANCE  =  ";F6;"  HZ" 

9216  FRIN/i  "DYNAMIC  ELECTFC^ECEAN ICAL  CCLPLING  COEFFICIE 
NT  =  " ; K 

^217  PRINT  "ELCCRED  RESISTANCE^  ";R;"CHMS" 

5231  PRINT  USING  92  32 

•?2  3  2  IMAGE  2/ 

9  250  FETURN 

9  30b  REM  **CALC-n/M** 

9  301  PRINTER  IS  70 1,76 

9  302  PRIN.T  @  PRINT 

9  30  3  PPINT  "vALLES   FOR  MEASUREMENTS  IN  WATER  FRCM  "  ;  3  $  ; 

"  DATA" 

9  304  PRINT 

S3 1 1  PRINT  "ELCCRED  INDUCTANCE  =  ";A£S(L0);"  HENRIES  AT 

";Fb;"EZ" 

9  312  PPINT  "ELECTPICAL  QUALITY  FACTOR  =  ";C1 

9  313  PRINT  "MECHANICAL  ^UALITY  FACTOR  =  ";£S 

9  314  PRIN.1  "ELECTPICAL  RESONANCE  =  ";F1;"  HZ" 

9315  PRIi'^T  "MECHANICAL  RESONANCE  =  ";F9;"  fiZ" 

931b  PRI:>.T  "DYNAMIC  ELECTROMECHANICAL  CCLPLING  COtFFICIH 

N  l  =  "  ;  K 

9  31c    PPIN.T    "POTENTIAL    EFFICIENCY    =    "  ;  P I 

9319    PRINT    "FREQUENCE    OF   OPTIMUM    EFFICIENCY    =■    ";,^;"     KZ" 

932U    PPINT    "MECHANICAL    fC^EZ    uTILI^AilCM    FACTOR" ;A£S (Rb ) 
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9322 
y323 
9  3  50 
9996 
9999 


PFIK'I    USING    9  32  3 

IftiAGE    3/ 

RETURN 

DISP    "THE    END" 

£NC 
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1  REfc  **PICKBW** 

2  OPTION  BASE  1 

3  SHORT  Q0,Q1,Q8,Q9 ,Q  ,E  ,E2  ,£ 3  ,E5 ,RS ,R7 ,X 2 ,X 3 ,X5 ,X6 ,X7 ,X8 
,G7,G8 ,E2 ,B3,B5,B6 ,hl ,E6 ,K,K1,P5,P6 ,a 

7  SHOFT  AfF2,F3,F4,F5,H3,H4,Jl,J2,J3,J4fG,B,R,X 

8  INTEGER  M,  X  ,  Yl ,  I ,  Y2  ,J  ,  Z2  ,Z  3  ,Hl  ,H2  fL  2  ,S  1 

10  DIM  B$[2]  ,A$[2]  rT$[l]  ,0(300,2)  ,Al$.[4]  ,D$[30C]  ,L$.(20]  , 
R$.[l]  ,I$[1]  ,S$.[10]  ,C$[5]  fG${23]  ,H$.[23] 

11  REAL  S  ,N,N8  ,N9  ,Rl  ,R9,X0,X9  ,L,C0  ,L0  ,0  ,P1,H8  ,H9  ,P,W0,G0 
,G9  ,E0 ,B9 ,N(10)  ,M(1Q)  ,F  0  ,  F  1  ,F6  ,F8  ,F  S 

14  FEAL  S8  ,E4,E6  ,L5fD  ,D1,D2  ,D3,D4 

15  X6,D1,D2  ,D3,D4  =  0 

16  Q,E2  ,E3,E5,RS ,R7 ,X2,X3,X5=0 

17  Z3,A,F2,F3,F4,F5  ,B6,£7,B&=0 

18  S8  ,E4  ,E6  ,M,Y,Y1  ,1  ,Y2,d  ,Z2=0 

19  K  ,K1,P5,P6  ,W,Q0,QI,C8,C.9  ,E  =  0 
2G  CO  ,L0  ,C,P1  ,B,Hti  ,H9  ,P,i\0  ,S=0 

21  N  ,N8,N9,R0  ,R9  ,X0  ,X9  ,GG  ,G9=0 

22  BO  ,E9  ,K1  ,H2  ,F0  fFl  ,F6  ,F8  ,F9=0 

24  DISP  "THIS  PROGRAM  IS  DESIGNED  FOF  AN  EXPERIENCED  CPE 
RATCR.  THE  BANDWID'IH  FOR  DATA  COLLECT." 

25  DISP  "AND  FREQ.  TO  GET  BLOCKED  DAI  A  MLL  BE  ThE  ChCIC 
E  OF  THE  OPERATOF.  "HYDRA2"'  PROGRAM  IS" 

26  DISP  "DESIGNED  TC  CCVEF  10  X  BW  AUTOMATICALLY.  HIT  "C 
ONT'  IF  READY  TO  PROCEED." 

27  PAUSE 

28  CLEAR  @  BEEP 

.25  DISP  "ENTER  TYPE  OF  TRANSDUCER  TC  BE  MEASURED?  ("ft  '  FC 
R  MAGNETIC  COUPLING  OR  "E"  FOR  ELEC.)" 

42  INPUT  T$£  DISP  " TC  GET  A  COMPLETE  SET  OF  DATA  YOU  NEE 
D  MEASUREMENTS  IN  BOTH  AIR  AND  WATER.  " 

45  DISP  "DC  AIR  FIRST. IN  WHAT  MEDIUM  ARE  YOU  OPERATING?  ( 
ENTER''!'  FOR  AIR  OR  '2'  FOR  WATER)" 

46  INPUT  M@  CLEAR  fc  BEEP 
4  7  IF  w=2  THEN  4  9 

40  CS="AIR"  y  GOTO  51 

49  C$="WAIER" 

51  PRINTEF  IS  2  ^  PLOTTER  IS  1 

55  DISP  "SET  THE  DRANEIZ  FOR  ADMI TTANCE  .SET  CN  LOWEST  SC 

ALE  THAI  WON'T  PEAK  DURING  EUU.lCCNT)" 

6C  PAUSE 

6  5  PRINT  $    PRINT  @  PRINT 

70  PRINT  "ADMTTANCE  IN  ";C$ 

75  GCSUE  7000 

80  IF  M=2  THEN  90 

85  F0  =  F6  (?  C0=Cj  &  GCTC  100 

90  F1=F6  (§  C1"C 

100  GCLEAR  @  CLEAR  (e  BEEP 

105  DISP  "SET  DRANETZ  FOP  IMPEDANCE- ( SCAL £  NEEDED)-HIT  C 

ONT  " 

in 


110  PAUSE 

115  PRINT  @  PRINT  @  PRINT 

120  PRINT  "IMPEDANCE  IN  ";C$ 

125  GOSUB  7000 

130  IF  M=2  THEN  140 

135  F8=F6  @  Q8=C.  @  GOTC  141 

140  F9=F6  (§  g9=G 

141  CLEAR    (?    BEEP 

142  PRINTER    IS    701,76    @    PRINI    USING    144 

14  4    IMAGE    2/ 

145  PRINT    "F(Y-AIR)=    "  ;  INT  (F  G)  ,  ,"F  ( Y-WAT)  =    "  ;  INI  (F  1)  ,  ,"  F 
(Z-AIR)  =    "  ;IN.T(F8)  ,  ,"  F  (&-aAT)  =  "  ;  INT  (F  9) 

146  PRINT    "£(Y-AIR)=    "  ;  <*  0  ,  ,"  Q  (  Y-WA1EP  )  =    ■  ;Ql ,  ,"Q  (Z-AIR)  * 
";Q8,,"Q(Z- WATER)"    " ;09 

147  PRINTER    IS    2 

148  CLEAR    @    DISP    "THESE    ARE    THE    INITIAL    ESTIMATES.  (HIT  "' 
CCNT'    WHEN    READY    TO    PROCEED) " 

149  PAUSE 

150  GOSUE  250  1 
152  CLEAR  @  EEEP 

15  3  DISP  "10  COLLECT  ADMITTANCE  DATA  (FOR  ELEC.  COUPLING 
)  ENTER  "1".  EN1ER  "2  '    FOR  IMPEDANCE" 

15  4  INPUT  L2@  PRINTER  IS  2 

155  IF  L2=2  THEN  158 

156  R$="G"  g  I$  =  "B"  is   S$="ADMITTANCE"  (<  H$="CONDUCTANCE 
(MICRCMHCS)"  £  G$="SUSCEPTANC£  (MICRCMHCS)" 

157  GOTC  159 

158  R$="R"     &    I$  =  "X"     fe    S$  =  "  IMPEDANCE"     I?    G$=" REACTANCE     (CH 
MS)"     §    H$="RESISTANC£     (OHMS)" 

159  IF    L2=l    THEN    176 

160  CLEAR  (§  DISP  "SET  ON  2  AND  ENTER  SCALE  FACTOR  TO  COL 
LECT  DATA . " 

161  INPUT  S&@  PRINI  S$;"  IN  ";C§  £  PRINT  "SCALE  =";£8;" 
OHMS" 

162  CLEAR 

164  GCSUD    8600 

165  IF    Y=2    THEN    152 
ioo    IF    M=2    THEN    16  9 

167    R8=E    (s    X3  =  £2    %    X8  =  E3    @    F4=£4 

16a    X6=E5    £    F5=E6    (5    Fb  =  El    @    GOTO    198 

169  R7  =  E    @    X2=E2    @    X7=E3    @    F4=E4 

170  X5=E5    @    F5=E6    (*    F9=E1    &    GOTO    198 

176  CLEAR    @    BEEP 

177  DISP    "SET    ON    Y    &    ENTER    SCALE    FACTOR    IN    MICRCMHCS    TC 
COLLECT    DATA." 

178  INPUI    S8(y    PRINI    S$;"    IN    ";C$    <s    PRINT    "SCALE     =";S6;" 
hlCROMHCS" 

179  GOSUE    8600 

160  IF    Y=2    THEN    152 

161  IF    M=2    THEN     185 
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162   G8=£    8    B3=E2    (e    £6=E3    @   F2=E4 

183    B6=E5    (9    F2=E6    8    FC=E1    @    GCTC    196 

165  G7=E    @    B2=E2    @    E7=E3 

166  F2=E4    @    £5=E5    @    F3=E6     8    F1=E1 
196   PRINT    USING    199 

199  IMAGE     2/ 

200  BEEP  @  DISP  "DG  YOU  DESIFE  A  LIST  OF  DA'IA  NEAP  RESON 
ANCE?  (IOC  PCINIS) (L=YES ,2=NC) " 

201  INPUT  Y 

2  02  IF  Y  =  2  THEN  300 

204  IF  L2=2  THEN  240 

210  PRIN.T  "ADMITTANCE  DA1A  FOR  "  ;C$;"  IN  MICROMHOS"  £  PR 

INT  USING  215 

215  IMAGE    2/  ,"  FREQUENCY"  ,  4X  /'REAL"  f6X  ^IMAGINARY" 

216  GOTO    265 

240  PRINT  "IMPEDANCE  DATA  FOR  ";C$;"  IN  OHMS"  @  PRIN1  US 

ING  215 

265  GCSUB  8100 

300  CLEAR  <s   BEEP  @  DISP  "DO  YOU  DESIRE  A  PLOT  OF  DATA?  (  1 

=  YES  ,2=NC) " 

30  5  INPUT  Y§  GCLEAR  £  CLEAR 

315  IF  Y«2  THEN  450 

320  DISP  "ENTER:  1=  G/R  VS.  FREw;  2=  B/X  VS.  FREO;  3=  £ 

VS.  G/X  VS.  R;  4=  END  PLOTTING  LOOP." 

325  INPUT  Y 

3  30  IF  Y  =  2  THEN  3  65 

335  IF  Y  =  3  THEN  425 

340  IF  Y  =  4  THEN  4  50 

345  GCSUE  6200 

356  LDIR  0  £  PEN  1  @  PENUP 

257  t^GVE  D(l,l)+.2*  (D(30Q,1)-D  (1,1)  )  ,-(E*.8) 

36  0  LABEL  R2>;"  VS.  FREQUENCY  IN  ";C$ 

361  GOTO  220 

265  GCSUE  8300 

2S7  LDIR  0  is  PEN  1  g  PENUP 

2£6  MOVE  D( 1,1)+ .2*  (D (30o,l)-D (  1,1) )  ,-L-.5*L 

4Cu  LABEL  1$;"  VS.  FREQUENCY  IN  ";C^» 

401  GCTC  320 

42  5  GCLEAR  g  CLEAR  £  GOSUB  6  4Ub 

431  FEN  1 

442  LDIR    C    (2    PENUP 

443  MOVE    -(. 5*L)  ,-L-.2*L 

445  LABEL  " INPU1  ELECTRICAL  ";S^;"  PLOT  ( " ; I  $ ; " vS. " ; E$; " 
)  IN  ";C$ 

446  GOTO  320 

450  IF  L2=2  THEN  645 

451  IF  T$="E"  THEN  525 

^5  2  GCLEAR  (§  CLEAR  @  DISP  "I  AM  CALCULATING  MCTICNAL  DAT 

A" 

465  C0*B*.Q0QG01/(2*PI*F6) 
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466  LO=l/( (2*FI*F6)~2*C0) 

470  FOR  1=1  TO  300 

475  D<I,2)=D(I  ,2).-G 

48  0  D(I  ,3)=D(I  ,3).-D(I,l)  *B/F6 

46  5  NEXT  I 

490  GOSUB  8700 

500  GOTO  54  2 

525  GCLEAR  @  CLEAR  <§  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 

A" 

532  C0=B*.O00001/(2*PI*F) 

535  FOR  1=1  TO  300 

536  D  (I  ,2)=B  (I  r2).-G 

537  D  (I,3)=DU  ,  3).-D(I,l)*B/F 

538  NEXT  I 

540  GOSUE  8700 

542  BEEP  @  DISP  "DO  YOU  DESIRE  A  LIST  OF  MOTIONAL  DATA  N 

EAR  RESONANCE? ( 1=YES ,  2=NO)" 

54  3  INPUT  Y 

545  IF  Y=2  THEN  6  04 

57  5  PRINT  "MOTIONAL  ADMITTANCE  DATA  IN  ";C?  £  PRINT  US IN 

G  560 

580    IMAGE     2/  ,"  FREQUENCY"  ,4X  ."REAL"  ,6X  .'"IMAGINARY" 

595    GCLEAR    <s  CLEAR 

6C0   GOSUE    8100 

604    CLEAR    @   BEEP    @    BEEP 

6C5    DISF    "WANT    A    PLOT? ( 1-tES  ,2*NO) " 

610    INPUT    Y2 

615  IF    Y2=2    THEN    755 

616  PLOTTER    IS    705    @    PEN    1 
620   GOSUB    8400 

628  LDIR    0    <e    PENUP    @    PEN    1 

629  MOVE    - (. 5*L)  ,-L-.2*L 

6  30  LABEL  "MOTIONAL  " ;S$;"  PLOT  FOR  ";C$ 
6  42  GOTO  7  55 

645  CLEAR  (?  GCLEAR  @  DISP  "I  AM  CALCULATING  MOTIONAL  DAT 
A" 

646  IF  T$="M"  THEN  65b 

647  FOR  1=1  TO  200 
b46  D( I  ,2) =E( I  ,2)-F 

649  C0=1/ (2*PI*F*X) 

650  LG=l/(  (2*PI*F6)  ~2*C0) 

651  D(I,3)=D(I  ,3).-2*PI*D(I  ,1)  *LC 
6  52  NEXT  I 

653  GOSUB  8700 

654  C0  =  A£S(C0)  (?  L0=AES(L0) 

655  BEEP  $    DISP  "DC  ¥CU  hANT  A  LIST  OF  THE  DATA  NEAR  RES 
ONANCE?( l=YES,2=NO)" 

656  INPUT  Y 

657  GOTO  670 

656  L0=X/ (2*PI*F6) 
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662  FOR  1=1  TO  300 

663  D  (I  ,2)=D(I  ,2).-R 

664  D(I,3)=D(I  ,3).-D  (I  ,  1)*2*PI*L0 

665  NEXT  I 

666  GOSUfi  8700 

66  7    BEEP    @    DISP    "    DO    YOU    WANT    A    LIST    OF    DATA    NEAP    RESONA 

NCE?(1=YES ,2=NG)" 

668    INPUT    Y 

670    IF    i=2    THEN    706 

66  5    PRINT    "MOTIONAL    IMPEDANCE    DATA    IN    " ;C$    &    PRINT    USING 

666 
666    IMAGE    2/," FREQUENCY" ,4X  , "REAL"  ,6X," IMAGINARY" 
700   GCLEAR    <t   CLEAR 

705  GOSUS    8100 

706  CLEAR    @    BEEP    @    BEEP 

710    DISP    "WANT    A    PLOT? ( 1=YES  , 2-NC ) " 
715    INPUT    Y2 

720  IF    Y2=2    THEN    755 

721  PLOTTER  IS  70  5  @  PEN  1 
725  GGSUB  8400 

737  LDIR  0  @  PENUP 

738  MOVE  -(.5*L)  ,-L-.2*L 

740  LABEL  "MOTIONAL  ";S$;"  PLOT  IN  " ; C$ 
75  5  GCLEAR  @  CLEAR 

756  PRINT  USING  7  57 

757  IMAGE  3/ 

800  DISP  "I  AM  DOING  CALCULATIONS  FOR  XOU.  PLEASE  BE    PAT 

IENT." 

80  3  IF  M=2  THEN  104  0 

810  IF  T$="M"  THEN  1000 

825  IF  L2=2  THEN  660 

640    C0=E1/A£S (E4-E6) 

045    C1=A£S(£2/(2*PI*F0) ) 

o55    GOTO    665 

obO    w8=El/A£S (E4-E6) 

665    C1=AES( 1/ (2*PI*F8*E2) ) 

6b5    K= 1- (FG/F6) ~  2 

855    K1=C  3/ (ABS(C0)+C1) 

S  30    GCSUB    9000 

931    PRINTER    IS    2 

94  0    GOTO    13  50 

10UU    K=l- (F6/F0) "2 
1005    IF    L2=2    THEN    1025 

1010    wU=El/ABS (E4-Eb) 

1015    GOTO    lu30 

1025    wo=El/AES(E4-E6) 

1l30  GCSUB  9200 

1031  PRIN1ER  IS  2  I?  GOTO  135u 

1040  IF  T$*"M"  THEN  1150 

lo4  1  BEEP  @  DISP  "IF  SOU  HAVE  LATA  FOR  AIR,  ENTEF  1;  ELS 
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E  2" 

10  42  INPUT  Y 

10  4  3  RAD 

LC44  IF  L2=l  THEN  1110 

1045  IF  Y=2  THEN  1062 

1060  P6=D4* (D3-D4)/(D3*R7) 

1061  WG=F3 

1062  £l=El/A£S(E4-£6)  @  K»1-(P  1/F9)  *2 
1070  GGSUB  9100 

10  76  PRINTER  IS  2  @  GOTC  13  50 

1110  IF  Y  =  2  THEN  1135 

1126  W0=F] 

1130  P6=D2* (Dl-D2)/(D1*G7) 

1135  Q1=E1/AE£(£4-E6) 

1140  K=l- (F1/F9) "2 

1145  GGSUB  9100 

114  6  PRINTER  IS  2  @  GOTG  135U 

1150  RAD 

1151  DISP    "IF    YCU    HAVE    DATA    FOR   AIR,    EN'lER     1;    ELSE    2" 

1152  INPUT    Y 

1153  IF  L2=l  THEN  1180 

1154  IF  Y=2    THEN  1162 

1155  fi=.5*ACS(  (R7-R9.)/D4) 

1156  H6=SIN(B)~2  @  H9.=COS(£)~2 

115  7    P5=(  (1+D3/R0*H9)  ~.5.-( 1-D3/RC*H8)  * . 5 )/ (  (  1+D3/R0*H 9)  ' 
.5+  (1-D3/RU*H8)  "  .5) 

1159  P=D4*SIN  (2*E)/  (4*R9*<^9)     fc    VvC=F9*  (  2*P  +  (P  '2+  1)  "  .5  ) 

1160  P6=D4* (D3-D4)/ (D3*R7) 

1162  »w9=£l/ABS(E4-E6)     @    K=  1-(F  9/F  1)  "2 

1163  PRINT    "B="  ;E,"K=";K  ,"  Kl  =  "  ;Kl,"H6="  ;rib  ,".49="  ;n9 

1164  GCSUB  9300 

1168  PRINTER  IS  2  y  GOTC  135  0 

1160  RAD 

1132  IF  Y  =  2    THEN  1210 

1164    E=.5*ACS( (G7-G9)/D2) 

1156    H6  =  SIN(B)~2    @    H9=CGS(B)~2 

1188    P5=(  (1+D1/G0*H9)  \ 5- ( L-D 1/C0* Ho )  \  5 )/ (  (  1+C 1/GU *H 9 ) 

.  5+  (1-Dl/GG*H8)  *  .5) 

11S0    P=D2*SIN  (2*B)/ (4*G9*<*1)     @    ;vO=Fl*  (  2*P  +  (P  ~2  +  1)  "  .  5  ) 

1195    Po=D2* (D1-D2)/ (D1*G7) 

1210  G;S=E1/A£S(E4-E6)     (e    a= L-(F 9/F 1)  ~2 

1211  PRINT    "£=" jBrMK*";KfM Kl*"?Kl ,MH8«" ;ho ,"H9=" ;H9 
1215    GCSUB    930o 

1225    PRINTER    IS    2    @    GOTO     1350 

1350    DISP    "DC    *CU    DESIRE   ANOTHER    IN    THIS    MEDIUM?  (.TO    GE'I 

THE    OTHER    TYPE    DATA?  )  (  1  =  YES  ,  2=i\C  )  " 

13  55    INPUT    Zl 

13  57  PLOTTER  IS  1 

13b0  IF  Zl  =  l  THEN  152 

136  1  PRINT  USING  136  2 
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13  62  IMAGE  3/ 

1365  DISP  "DC  YCU  DESIRE  A  RUN  IN  ANOTHER  MEDIUM? ( 1= YES , 

2*80)" 

137U  INPUT  Zl 

1375  IF  Zl  =  l  THEN  45 

1360  DISP  "JUST  FOR  THE  RECORD,  INPUT  WATER  TEMP,  AIR  TE 
MP,  TRANSDUCER  SER  NO,  MODEL  NO.  ," 

1361  DISP  "AND  TYPE  (E  OR  M)" 

138  5  DISP  "INPUT  "'0"  IF  INFO  IS  uNKNGWN" 
13  90  INPUT  W,A,S  ,N  ,L$ 

13  51  PRINTER  IS  701,76 

1400  PRINT  "WATER  TEMP  =  "  ;  W  ,  ,"  AI  P.  TEMP="  ;A  ,  ,"  SER.  NO.  =  ";S 
,  ."MODEL  NO.  =  "  ;N  ,  ,"IYPE  ";L$ 

14  10   GOTO    9996 

2501  REM    *SHUNT* 

2502  CLEAR    <5    DISP    "INPUT   A    PERCENTAGE    OF    RESCNANT    FREQ. 
TO    USE   TO    GET    SHUNT    (BLOCKED)     VALUES." 

2  50  3    DISP    "INPUT    A    VALUE    FROM    0    TO    1." 

2  50  4    INPUT    Bl 

2  513  CLEAR  £  BEEP  @  DISP  "vvE  ARE  FINDING  SHUNT  VALUES  FO 

R  G/B  OR  R/X." 

2515  A$="FR"  @  B$  =  "HZ"  @  S 1=  1 

2519  IF  T$="M"  THEN  2525 

2520  IF  M=2  THEN  2522 

2521  F=F0*B1  @  GOTO  2535 
2  5  22  F=F]*B1  is    GOTO  25  3  5 

2525  IF  M=2  THEN  2527 

2526  F=F8.-S1*F6  g  GOTO  2535 

2527  F=F9-B1*F9 

2535  D$=VAL$  IF.) 

2537  OUTPUT  717  ;A$,D$,£$ 

2536  WAIT  1000 

2540  DISP  "SET  DRANETZ  £RE^.    SCALE  TO  CCvEF  FREw .  ON  SYN 

THESIZER  .ZERO  METERS  .SET  ON  FS.(CONT)" 

254  1  PAUSE 

254  2    DISP    "HERE    yvE    GET    THE    NORMALIZATION    FACTORS." 

2543  OUTPUT  709  ?"VT1" 

2548  GCSUB  8549 

2550  G  ,B  ,R,X=0 

2551  A$="FR"     £    Ei?  =  "HZ"    &    D$»VAL$(F) 
2553    OUTPUT    717    ;A$,D$,E$ 

2561  CLEAR    @    DISP    "SET    DRANETZ    ON    i  ,SET    SCALE, SET    CN    i^CR 
M    FOR    FILTER    +    PHASE, PLUG    IN    TRANSDUCER." 

2562  EEEP    @    DISP    "ENTER    SCALE    FACTOR    1l\    MICRGKHCS.SAME    S 
CALE   AS    FOR    DATA. (WE    GET    SHUNT    G/B)" 

2564  INPUT    S£(?    CLEAR    g    DISP    "I    AM    WORKING    TO    GET    SHUN'I    v 
ALuES"     £    SEcP 

2565  FOR    1=1    TO    10 

2  5  66    OUTPUT    717    ;A:?,D$,B$ 
2567    rtAIl    1000 
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25bb  OUTPUT  709  ;"AC3" 
2  56  9  ENTER  7  0S  ;  M(I) 

2570  G=G+h(I) 

2571  OUTPUT:  709  ;"AC4" 

2572  ENTER  70S  ;  N(I) 

2573  B=B+N(I) 

2574  NEXT  I 

2575  G=G/(I-1)  @  E=E/(I-J) 

2576  G«G*S8/N8  g  8*8*S8/N9  @  BEEP 

2  576  DISP  "SET  CRANETZ  ON  o ;  SET  SCALE  FOR  MAX  RESPONSE.  N 

OVE  TRANS.  INPUT.  ENTER  SCALE  FACTOR" 

257S  INPUT  Sot  CL££E  @  DISP  "I  Kb    WORKING  TO  GET  SEUNT  V 

ALOES" 

2580  EOF  1=1  TO  L0 

2  561  OUTPUT  717  ;A$,D$,5$ 

2562  aAIT  1000 

258  3  OUTPUT  709  ;"AC3" 

2584  ENTER  70S  ;  K(I) 

258  5  R=R+M(L) 

2566  OUTPUT  709  ;"AC4" 

2  5b' 7  ENTER  7uS  ;    N(I) 

2568  X=N(I)+X 

256  9  NEXT  I 

2590  R=R/(I-1)  @  X=X/(I-J) 

2591  R=R*S6/N8  @  X=X*S6/N9 

2592  IF  T$«"M"  THEN  2610 
2595  IF  N=l  THEN  2600 

2  5  96  R9=R  £  X9=X  fc  C9=G  @  3S=£  @  GOTO  2601 

2600  G0=G  £  LU=B  (g  R0=R  fe  X0=X 

2601  PRINT  "SHUNT  VALUES" 

2602  PRINT  USING  2603 
2o03  LMAGE  1/ 

260  5  PRIN1  "G0  =  ";G0  ,  l"E0  =  ";BC  ,  #"R0=";RU  ,  ,"X0  =  ";a0  ,  ,"GS=" 
;GS,  ,"£9  =  ";  £9,  ,"RS=";R£I  ,"X9=";X9 
2606  PRINT  £  GOTO  26  20 

2610  IF  Sl=2  THEN  262u 

2611  J1=G  t  u2=£  @  J3  =  R  g  J4=X 

2612  IF    M*2    THEN    2b 15 

2613  F.=F8  +  E1*F&    @    Sl  =  2    @    GOTO    255u 
2615    F=F9+B1*F9    @    Sl=2    fe    GOTO    255o 

26  20    G=.5*(G+J1)    (e    o=.5*(B+J2)     t   R*.5*(R+J3)     @    X*.5*(J4+ 

X)     %    GOTO    25  9  5 

2  6  30  GCLEAR  @  CL££R  £  BEEF 

2635  CISP  "ENTER  '  1* '  IF  ALL  IS  WELL;  ''2'  IF  YOU  NEED  A  B 

EPEAT." 

2b 40  INPUT  Y 

2645    IF    Y=2    THEN    250  1 

2  6  50    RETURN 

7000    REN    **F-<*** 

700  1    GCLE£P     g    CL££R 
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7005  DISP  "ENTER  LOWER  ,  uPPER  FREC.  FOR  SWEEP. IF  RESCNAN 
C£  UNKNOWN  USE  10000,2  00000.  (FlfF2)" 

7010  INPUT    F4,F2 

7011  DISP    "ENTER    RMS    VOLTAGE    NEEDED.  (30  00    MV    FOR    LOW    FRE 
U_.    AS    LOW    AS    2C0C    MV    FOR    HIGH)" 

7012  INPUT   Ale   Al$  =  VAL$(A) 
7020    CLEAR    @    BEEP 

7025    A$=HAM"    $    B$="MR" 

70  35  OUTPUT  717  ;A$,A1$,B$ 

7040  OUTPUT  709  ;"AC3VT3" 

7045  F3=CEIL(  (F2-F4)/30C) 

7  06  0   CLEAR    %    CISP    "I    AM    WORKING    TO    GET    THE    SPECTRUM    DAI A 

FOR   YOU, .THEN    WILL    SHOW    &    MAKE    A    PLOT" 
7065    FOR    1=1    TO    300 
7  0  70    D  (I  ,L)=F4  +  I*F3 
7075   A$="HZ"    §   B$="FR" 
7083    D$=VAL$(D (I ,1) ) 
7005   OUTPUT    717     ;B$,D$,A$ 
7090    OUTPUT    709     ;"AC3VT3" 
7095    WAIT    100 
7100    ENTER    709    ;    D(I,2) 
71u5    NEXT    I 

7110    F5=F4-.1*(F2-F4)     @    F6=F2+ . 1* (F2-F 4) 
7115    F6=F2-h.l*(F2-F4) 
7119    PLOTTEF    IS    1    €    GCLEAR    §    CLEAR 
7125    SCALE    F5  ,F6 ,- . 1  ,1 . 2 
7130    XAXIS    0,200C,F4,F2 
7135    tfAXIS    F4,  .1  ,-.1  ,1.2 
714  0    FOR    1=1    TO     30  0 
7145    PENUP 

7150    PLOT    D(I  ,1)  ,D(I  ,2) 
7155    NEXT    I 
7I5c    GR£PH    f>    CCtY 

7  160    DISP    "NEED    ANOTHER    RUN?     (  1=1E£  ,  2  =  NC ) " 

7161  INPUT    Y 

7162  IF    Y*l    THEN    7005 

716  3    CLEAR    &    GCLEAR    @    BEEP 

7165    DISP    "ENTER    DECISION    POINT    FOR    AMPLITUEE     (0    TO    1.2) 

.  (LESS    THAN    THE    MAX    DISPLAYED) 

7170    INPUT    S£   CLEAR    fe    BEEP 

"/1C0    PRINT     "AMPLITUDE    IN    vCLTS";  "         FREQUENCY" 

7161    PRINT 

7  190    FOR    1=1    rC    30  0 

7195    IF    D(I,2)<S    THEN    7210 

72  00    PRINT    USING    7  2U5    ;    D ( I , 2) ,D ( I , 1) 

7205    IMAGE     IX  ,D.DDCDCDDC , 1CX  ,DDODUC .DD 

7  2 lu    NEXT    I 

7  215  PRINT  USING  7  22o 

72 2U  IMAGE  3/ 

7225  CLE^R  it  EEEP 
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7226  DISP  "IF  YOU  ARE  READY    TC  CONTINUE  ENTEF  "l";  OTHER 
WISE  ENTER  "1 ' . " 

7227  INPUT  Y 

7228  IF  Y=2  THEN  7160 

7255  DISP  "WHAT  IS  CENTER  FREQ  ,  HALFWID1H  TC  CONSIDER?" 

7236  INPUT  F,N(?  CLEAR  £  BEEP 

72  40  DISP  "I  AM  WORKING  TO  GET  GOOD  FRE<*  AND  Q    DATA" 

7241  FOR  1  =  1  TO  30G 

7242  D  (I  ,  1)  =F.-N+I*N/150 

7243  A$="H2"  @  B$="FR" 
7245  D$=VAL$(D (I ,1)  ) 
7255  OUTPUT  717  ;£$,D$,A$ 
726C  OUTPUT  709  ;"AC3VT3" 
7265  wAIT  100 

7270  ENTER  7o9  ;  D(I,2) 

727  5  NEXT  I 

7  28  5  CLEAR  @  GCLE£R  i§  DISP  "I  AM  FINDING  THE  ABSOLUTE  MA 

X  AND  FREQ  UPPER  AND  LOWER" 

7290  E1,E4=C  @  HfHl,H2,H3,H4=l 

7295  62,E3=50 

7300  FOR  1=2  TO  300 

7310  IF  D(I ,2)  <D(H,2)  THEN  7  340 

7315  IF  D(I  ,2) =  D(H,2)  THEN  7325 

7325  A6=D(I,2)  £  F6=C(I,1)  (9  H=I 

7340  NEXT  I 

7345  A7=A6/SQR(2) 

7350  FOR  1=1  TO  K 

7355  IF  A7=D(If2)  THEN  7405 

73b0  IF  A7<D(I,2)  THEN  7385 

7  36  5  IF  D(I,2)<B]  THEN  7425 

7370  Bl=D(I,2)  @  H1=I  @  GOTO  7425 
7385  IF  D(I,2)>82  THEN  7400 

7390  B2=D(I  ,2)  @  H2  =  I 

7  40  0  GOTO  7  4  25 

7405  B1,B2=D(I,2)  @  Hi  ,H2=I  @  F7  =  D(H2,1)  @  GOTO  7440 

742  5  NEXT  I 

7430  X=(A7-E1)/  (B2-E1) 

7  4  35  F7  =  X*  (D(d2  ,  l).-B(Hl,l)  )+D  (ril,l) 

7  4^0  FOR  I=H  TO  30 C 

7445  IF  A7=D(I,2)  THEN  7495 

7  4  50  IF  A7>D(I,2)  THEN  7^75 

7  455  IF  D(I,2)>B3  THEN  7470 

7460  E3  =  D(I  ,2)  fe  H3  =  I 

7  4  70  GOTO  7  515 

7475  IF  D(I,2)<=B4  THEN  7490 

7480  B4  =  D  (1,2)  g  H4  =  I 

7490  GOTO  7515 

7495  B3,E4=D(I,2)  @  ri3,H4=I 

75u5  L5=C  (H3,  1) 

7  510  GOTO  7  5  30 

120 


7515  NEXT  I 

7520  X=(A7-E4)/(E3-E4) 

7  525  L5=-(X*  (D(H4  ,1).-D(H2,1)  )  )+D  (H4.1) 

7530  Q=F6/(L5-F7) 

7535  PRINT  "CENTER  FREQ  IS  ";F6 

7  5  36  PRINT  "^    IS  " ;  Q 

7538  DISP  "IF  YOU  APE  READY  TO  PROCEED,  ENTER  "l'.  TO  RE 

RUN  FOR  BETTER  VALUES  ,  ENTER  '2'." 

7  53  9  INPUT  Y 

7540  IF  Y  =  2  THEN  7235 

75  45  RETURN 

6100  REM  **DATA  LIST** 

8101  PRINTER  IS  2 

8102  FOR  1=99  TO  199 

810  3    PRINT    USING    8  10  5    ;     D  (I  ,  1)   ,D  ( I  ,2)  ,D  ( I  ,  3) 
81U5    IMAGE    DDDDDD  .D  ,2X  .D.DDDE       ,2X, D.DDDE 

8106  NEXT  I 

8107  PRINT  USING  8106 
6108  IMAGE  3/ 

6110  RETURN 

6200  REM  **PLOT  RE** 

8201  GCLEAR  @  CLEAR 

6  202  H8=D  (1,1) -.15*  (D(  300,1)  -D(  1,1)  ) 

620  3  H9  =  D (3CG,l)  +  .l*  (D  (  300,1).-D(  1,1)  ) 

6205  PLOTTER  IS  7U5  &  PEN  1 

8209  SCALE  H8  ,H9  ,- ( . 8*£ )  ,E  +  . 25*E 

6210  XAXIS    0 , (D( 300 ,1)-D (1,1) )/10  ,D  (  ]  ,])  ,D( 300  ,1) 

6211  YAXIS    D( 1,1)  ,E/5  ,-(. L*E)  ,£+.2*E 
8  212    PEN UP 

8  215    FOR    1=1    TO    30  0 

6216    PLOT    D(I ,1)  ,D(I  ,2) 

8217   NEXT    I 

6216    PENUP    (y    DEG    @    LDIR    C,SIN(9b) 

0221    FOR    L1  =  D(1,1)    TO    D(2uG,l)     STEP    (D  (  20 0  , 1 ) -C  ( 1 ,  1)  ) /  10 

@     PENUP 
8  222   MOVE    LI,-  (.  16*E) 
8  22  2    LAbEL    INT (LI) 
6  224    NEXT    Ll 
o226    LDIR    0    @    PENUP 

8227    I^.CvE    D(l  ,1)  +  .2*  »D  (  300  ,1)-D  (1  ,1)  )  ,-(  .7*E) 
6226    LAEEL    "FREQUENCY     (HZ)" 
8230    PENuP    £    PEN     1 
d232    LDIR    0 

b234    FOR    LL=-(.8*E)     TO    E+.2*E    STEP    £/ 5    fe    PENUP 
62  3  5    MOVE    D  (1,1) -.09*  (D  (  2  00  , 1  ).-D  (1,1))  ,L  1 
6  2  36    LAEEL    INl(Ll) 
b2  37    NEXT    Ll 
624b    LDIR    0,SIN(9Q) 

6241  MOVE    D(l,l).-.  1*  (D  (3Co,l)-D  (1,1)  )  ,.1*E 
o24  3    LAEEL    fi§ 
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8 2 50    RETURN 

8  200    REM    **PLOT    IM** 

8301   GCLEAR    @   CLEAR 

8  30  3   H8  =  D(  1  ,1).-.  15*  (D  (300,1  )-D  (  1,1)  ) 

8  30  4    H9=D(30  0  , 1)  +  .1*  (D  (  300  ,  L)  -0  (  1 , 1)  ) 

8305   L=MAX(ABS(E5.)  ,ABS(E3)  ) 

8  30  6  PLOTTER  IS  70  5  @  PEN  1 

8312  SCALE  H8  ,H9  ,-L-. 5*L  ,L+ . 25*L 

8313  XAXIS    0,(D(  300,1)-D(1,1)  )/10,C  (  1  ,1)  ,D(.30C  ,1) 
6314    YAXIS    D(  1,1)  ,L/5  ,-L-.2*L ,L+.2*L 

8315   PENUP 

8  318   FOR    1=1    TO    20  0 

8  325   PLOT    D(I,1)  ,B(I,3) 

8330  NEXT    I 

8331  FENUP    @    DEG    §    LCIR    0,SIN(90) 

8333    FOR    L1  =  C(1,1)    TO   D(.200,l)    STEP    (D  (  200  , 1 )  -D  ( 1  ,  1)  )  /  10 

@    PENUP 
8225   MOVE    L1,-L-.2*L 
6  3  36    LABEL    INT (LI) 
8337   NEXT    LI 
8340    LOIR    0    (?    PENUP 

8  341   MOVE    D(l,l)  +  .2*  (D(300,l)-0  (  1,1)  )  ,-L-.4*L 
6342    LABEL    "FREQUENCY    (H£)" 
8345    PENUP    @    PEN     1    <?    LDIR    0 

8352    FOR    LL=-L-.2*L    TO   L-h.2*L    STEP   L/5    &    PENUP 
6  3  55    MOVE    D  ( 1 , 1) -.  09*  (D  (.200  , 1) -D  (  1  ,1)  )  ,L1 
6  3  56    LABEL    INT(Ll) 
8  257   NEXT    LI 

6260    PENUP    (?    PEN     1    @    LCIR    C,SIN(90) 
8261   MOVE    D(l,l)-. 1*  (D (300,1)-D (1,1) )  ,-(L*.5) 
8  266    LAEEL    G$ 
8  260    RETURN 
8400    REM    **CIRCLES** 
6401    GCLEAR    @    CLEAR 
8402    L1=MAX(AES (E3)  ,AES(E5) ) 

6403  L=MAX(A6S(E)  ,L1) 

6404  PLOTTER  IS  70  5  ig  PEN  1 

8406  SCALE  1.  25*  (-L-  .  1*L  )  ,  1 .  25*  (L+  .  1*L  )  ,-L-  .  2*L  ,L-h  .  15*L 

8407  XAXIS  0  ,L/5  ,-L-. 1*L  ,L+ .  1*L 
6406  YAXIS  0  ,L/5  ,-L-.  1*L  ,L+.  1*L 
8  409  PENUP  £  PEN  1 

34  11  FOR  1=1  TO  20  0 

6  412  PLOT  L  (I  ,2.)  ,0(1,2) 

6413  NEXT  I 

6414  PENUP 

84  15  LDIR  0  (9  PEN  1  0  PENUP 
8416  MOVE  -(L*.5)  ,-L-.l*L 
641o  L£EEL  H$ 
6420  PENUP  @  PEN  1  @  CEG 
d425  LDIF  0,SIN(90) 
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6  426  FCR  LL=-L-.1*L  TC  L+.1*L  STEP  L/ 5  @  PENUP 

6427  toCVE  Ll,-L 

6428  LABEL  INT (LI) 
6  429  NEXT  LI 

8435  LDIR  G  @  PEN  1  @  PENUP 

6436  FOR  L1=-L-.1*L  TO  L+.1*L  STEP  L/5  @  PENUP 

8437  MOVE  - (.866*L)  ,L  1 

8438  LABEL  IN.T(Ll) 
843  9  NEXT  LI 

8445  LDIR  G,SIN(9G) 

8446  fcOVE  -  ( .  95*L)  ,- ( . 3*L) 
6448  LABEL  G$ 

6450  RETURN 

8  5*9  REM  **NCRM** 

6550  N8  #N9=0 

6551  FOR  1=1  TC  10 
6  557  WAIT  500 

8556  OUTPUT  709  ;"AC3" 
8559  ENTER  709  ;  N(I) 
6560  N8=N8+N(L) 
6565  OUTPUT  709  ;"AC4" 
6  566  ENTER  70  9  ;  M(I) 

6567  N9=N9  +  fc(l) 

6568  NEXT    I 

8569  N8=N8/(I-1) 

8570  N9=-(N9/(I-1) ) 

8  57  2   PRINT    USING    8  573 

6  57  3    IiviAGE    3/ 

6575    RETURN 

6600    REN    **vALUES** 

66  01   CLEAR    9    BEEP 

66 05  IF    L2=l   AND    M=2    THEN    66  10 

6606  IF    L2=l   AND    K»l    THEN    6611 
66  07    IF    L2=2    AND    M=l    THEN    6b 12 
o6  06    IF    L2=2    AND    fo*2    THEN    66  12 
66  10    F6=F1    @    Q*Q1    @    GOTO    8615 
8611    F6=F0    @    C=w0    @   GOTO    66  15 
86  12    F6=F6    @    Q=^b    @    GOTO    do  15 
66  13    F6=B9    G   v=v9 

ob  15    DISP    "PLEASE    SELECT    THE    HALFWIDTH    OS    INTEREST    TC    CO 

LLECT    DATA.(ENTEE    HALFWIDTH) " 

oblb    INPUT    inj^    CLEAR    @    DISP     "I    AM   COLLECTING    THE    DATA." 

6625    FCR    1=1    TO    300 

b629    D (I ,1)=F6-N+I*N/150 

b630   A$="FR"    i   B$*"HZ" 

6632    D$=vAL3(D( I  ,1)  ) 

663;    CLTFLT    717    ;A$,D$,E$ 

6634    WAIT    250 

ob3  5   OOTPul    709     ;"AC3vT3" 

6636    ENTER    709    ;     Ho 
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8637  D(I  ,2) *H8*S8/N8 

66  38  OUTPUT  709  ;"AC4VT3" 
66  3  5  ENTER  70S  ;  H9 

8640  D  (I  ,3)=H9*S8/N9 

8641  NEXT    I 

8662   CLEAR    (§    DISP    "I    AW    FINDING    MIN/MAX    VALUES    &    ASSOC. 
FREQS" 

6663  GCSUE    8700 

6664  IF    L2=2    THEN    8687 

8665    PRINT    "GM£X  =  "  ;E  ,  ,"  FGMX=";E1  ,  ,"EGMX="  ;E2  ,  ,"EMX  =  "  ;£3  , 

,"  FBMX*"  ;  E  4  ,  ,"  BMI«n  ;E5  ,  ,"  FBMI=";E6 

56  8  6   GCTC    8  6  97 

6667    PRINT    "RKAX="  ;E  ,  ,"  FRMX=";E1  ,  ,"XRKX="  ;E2  ,  ,"  XMX  =  "  ;E3  , 

,"FXMX  =  "  ;E4,  ,"XMI=";E5  ,  ,"FXKI=";E6 

8697  DISP  "If  ALL  IS  kELL  ,  ENTER"!'.  ENTER  "2"'  TO  RETAK 
E  DATA." 

8698  INPUT  ¥ 
RET LBN 

REM  **MAXMIN** 
H,H1,H2=50 

E=D(50,2)     9    Ei,E4  ,E6=D(50  ,1)     @    £2  ,E  3  ,E  5=D  |5C  ,  3) 
FOR    1=51    TC    250 
IF    D(I  ,2)  <D  |Hr2)    THEN    6725 
IF    D(I ,2)=D (H,2)    THEN    8  72  0 
£=D(I,2)     £    £1=0(1,1)     tf    £2=C(I,3)     g   H=I 
IF    D  (I ,3) <D  (HI  ,3)    THEN    8740 
IF    D  (I , 3) =D  (HI  ,3)    THEN    67  35 
E3=D(I,3)     &    E 4=0(1,1)     (£   fl  1=1 
IF    D  (I  ,3)  >  =  D(H2,  3)    THEN    8751 
H2=I    <e   E5  =  D(I,3)     @    E6.=D(I,1) 
NEXT    I 
J1=D(H,2)     §    J2=(D  (H+1,2)-D(H-1  ,2)  )/2    @    J:=(D(H  +  1,2) 

+D(H-l,2).-2*JL)/2 

F6=-(J2/(2*J3) ) 

A=J  1+J2+F8+J  3*F6~2 

F6=D  (H.D+F6*  (D(H  ,1)-D  (H-l  ,1)  ) 

El=F6    §    e=A 

D.=A£S(£3-E5) 

IF    t-'=2    THEN    67  9  3 

IF    L2=l    THEN    6792 

03=0    @    GCTC    o7£6 

D1  =  D    is    GCTC    67S6 

IF    L2=l    THEN    8  79  5 

D4=L    @    GCTC    6796 

02=0 

PRINT    "Dl="  ;D1,  ,"D2=";C2  ,  ,"C3  =  ";D3  ,  ,"  C4«H;04 

RETURN 

REfc    **CALC-A** 

PRINTER    IS    70  1 ,76 

PRINT    ilj    PR  INI 
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6699 

6700 

8703 

8705 

8706 

6710 

8711 

8720 

8725 

8726 

6  7  35 

6  7  40 

8745 

8751 

8760 

+D(H- 

8  752 

8  75  5 

6770 

8775 

&786 

8  7o  9 

6790 

6791 

8  792 

b793 

8794 

8795 

67b»6 

8799 

90U0 

3uOi 

you2 

9003 

PRINT 

9004 

PRINT 

SOU 

PRINT 

;F;"HZ" 

9012 

PRINT 

9013 

PRINT 

9014 

PRINT 

9015 

PRINT 

9016 

PRINT 

NT  = 

B;K 

9017 

PRINT 

9021 

PRINT 

9031 

PRINT 

9032 

IMAGE 

9  0  50 

RE1URN 

9  IOC 

REM  ** 

9101 

PRINTE 

9102 

PRINT 

9103 

PRINT 

"  DAI A" 

9104 

PRINT 

Sill 

PRINT 

'HZ" 

9112 

PRINT 

9113 

PRINT 

9114 

PRINT 

9115 

PRINT 

9116 

PRINT 

MT  = 

";K 

9119 

PRINT 

"VALUES  MEASURED  IN  AIR  FROM  ";S$;"  DATA" 

"ELOCKED  CAPACITANCE*  " ;AES(C0);"  FARADS  AT  " 

"ELECTRICAL  QUALITY  FACTOR  =  ";Q0 
"MECHANICAL  QUALITY  FACTOR  =  ";w6 
"ELECTRICAL  RESONANCE  =";F0;"  HZ" 
"MECHANICAL  RESONANCE  =  ";F8;"  HZ" 
"DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 

"STATIC  COUPLING  COEFFICIENT*  ";K1 
"BLOCKED  RESISTANCE*" ;R0; "OHMS" 
USING  9  0  32 
2/ 

CALC-w** 

R  IS  701,76 

@  PRINT 

"VALUES   FOR  MEASUREMENTS  IN  WATER  FRCM  ";S$; 

"ELOCKED  CAPACITANCE  =  ";AES(CG);"  FARADS  AT 

"ELECTRICAL  QUALITY  FACTOR  =  ";CJ 
"MECHANICAL  QUALITY  FACTOR  =  ";QS 
"ELECTRICAL  RESONANCE  =  ";F1;"  HZ" 
"MECHANICAL  RESONANCE  =  ";F9;"  HZ" 
"DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 

"FREQUENCY  OF  OPTIMUM  EFFICIENCY  =  ";v\0;"  HZ" 


9120  PRINT  "MECHANICAL  POWER  UTILIZATION  FACTOF" ; AES (P6 ) 


9121 

PRINT  "  STATIC  COUPLING  COEFFICIENT*" ; Kl 

9122 

PRINT  "ELOCKED  RESISTANCE  =  " ; R9  ;  "OHMS " 

9132 

PRINT  USING  9  13  3 

9133 

IMAGE  3/ 

9150 

RETURN 

920U 

REM  * *CA LC -A /M** 

9201 

PRINT  EF  IS  7  0  1,76 

9  202 

PRINT  £  PRINT 

9203 

PRINT  "\/ALUES  MEASURED  IN  AIR  FROM  ";SS>;"  DATA" 

9204 

PR  IN 'I 

9211 

PRINT  "ELOCKED  INDUCTANCE*  "  ; AcS (LC )  ; "  HENRIES  A'l  " 

;Fb;' 

'HZ" 

y  212 

PRINT  "ELECTRICAL  QUALITY  FACTOR  =  " ; C;0 

9213 

PRIN1  "MECHANICAL  QUALITY  FACTOR  =  ";wb 

9214 

PRINT  "ELECTRICAL  RESONANCE  =";F0;"  HZ" 

9215 

PRIN'i  "MECHANICAL  RESONANCE  =  ";Fb;"  HZ" 

92  lt> 

PRINT  "DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 
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Nl    =    "  ;  K 

5217    PRINT    "ELCCKED    RESISTANCE*    " ; F 

92  31   PRINT    USING    92  22 

9222    IMAGE    2/ 

9  250    RETURN 

5  30  0    REM    **CALC-W/M** 

9 301    PRINTER    IS    7  01,76 

9  30  2    PRINT    @    PRINT 

9303    PRINT    "VALUES       FOR   MEASUREMENTS    IN    K.A1ER    FRCM    ";S$; 

"    DATA" 

9  304    PRIN'I 

5311    PRIN'I    "LLGCKED    INDUCTANCE    =    "  ;A£S(Lu);"    HENRIES    AT 

";F6;"HZ" 

9  312  PRINT  "ELECTRICAL  QUALITY  FACTOR  =  a;Ql 

*312  PRINT  "MECHANICAL  QUALITY  FACTOR  =  ";VS 

9314  PRINT  "ELECTRICAL  RESONANCE  =  ";F1;"  HZ" 

5315  PRINT  "MECHANICAL  RESONANCE  =  ";FS;"  HZ" 

9316  PRINT  "DYNAMIC  ELECTROMECHANICAL  COUPLING  COEFFICIE 

NT  =  " ; K 

9  31b    PRINT    "POTENTIAL    EFFICIENCY    =    " ;P5 

931b>    PRINT    "FREQUENCY    OF    OPTIMUM    EFFICIENCY    =    ";v*0;"    HZ" 

y320    PRINT    "MECHANICAL    PC^EP    UTILIZATION    FACTOR" ;AES (P6) 

5322    PRINT    USING    9322 
9  32  3    IMAGE    3/ 

5  3  50    RETURN 

S55b    DISP    ".THE    END" 

5555    END 
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